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Abstract 

We consider dynamic algorithms for maintaining Single-Source Reachability (SSR) 
and approximate Single-Source Shortest Paths (SSSP) on n-node m-edge directed graphs 
under edge deletions (decremental algorithms). The previous fastest algorithm for SSR 
and SSSP goes back three decades to Even and Shiloach [JACM 1981]; it has 0(1) query 
time and 0{mn) total update time (i.e., linear amortized update time if all edges are 
deleted). This algorithm serves as a building block for several other dynamic algorithms. 
The question whether its total update time can be improved is a major, long standing, 
open problem. 

In this paper, we answer this question affirmatively. We obtain a randomized algorithm 
with a total update time of 0 (min(TO^/®n^/^+°^^\ = O(mn®no+o(i)^ 

SSR and (1 -I- e)-approximate SSSP if the edge weights are integers from 1 to W < 2^°® 
and e > 1/ log° n for some constant c. We also extend our algorithm to achieve roughly the 
same running time for Strongly Connected Components (SCC), improving the algorithm 
of Roditty and Zwick [FOCS 2002]. Our algorithm is most efficient for sparse and dense 
graphs. When m = 0(n) its running time is and when m = 0(n^) its 

running time is O(n^+^n-i-o(i) Por SSR we also obtain an algorithm that is faster for 
dense graphs and has a total update time of 0 ( 771 ^/^ 71 "*/^+°^) _|_ 77 ^ 3 / 7 J.Jl 2 / 7 -l-o(l)^ which is 
0(n^+^/^) when m = 0(n^). All our algorithms have constant query time in the worst 
case and are correct with high probability against an oblivious adversary. 
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1 Introduction 


Dynamic graph algorithms are data structures that maintain a property of a dynamically 
changing graph, supporting both update and query operations on the graph. In undirected 
graphs fundamental properties such as the connected, 2-edge connected, and 2-vertex 
connected components as well as a minimum spanning forest can be maintained very quickly, 
i.e., in polylogarithmic time per operation [HK99, HLTOl, ThoOO, KKM13], where an 
operation is either an edge insertion, an edge deletion, or a query. Some of these properties, 
such as connectivity, can even be maintained in polylogarithmic worst-case time. More 
general problems, such as maintaining distances, also admit sublinear amortized time per 
operation as long as only edge deletions are allowed [HKN14a, HKN14b]. 

These problems when considered on directed graphs, however, become much harder. 
Consider, for example, a counterpart of the connectivity problem where we want to know 
whether there is a directed path from a node u to a node v, i.e., whether u can reach v. In 
fact, consider a very special case where we want to maintain whether a fixed node s can 
reach any node v under edge deletions only. This problem is called single-source reachability 
(SSR) in the decremental setting. It is one of the simplest, oldest, yet most useful dynamic 
graph problems. It was used as a subroutine for solving other dynamic graph problems 
such as transitive closure [HK95, Kin99, BHS07, RZ08] and strongly connected components 
(SCC) [RZ08]. Furthermore it is a special case of decremental (approximate) single-source 
shortest paths (SSSP) [ES81, BRll, HKN14b], which in turn was used as a subroutine for 
solving dynamic all-pairs shortest paths [Kin99, BHS07, BHS03, RZ08, RZll, Berl3, HKN13, 
ACT14]. Yet, no algorithm with sublinear update time was known for this problem. This 
was also the case for decremental s-t reachability (stR), where we want to maintain whether 
the node s can reach the node t, and s-t shortest path (stSP), where we want to maintain 
the distance from s to t. 

Related Work. To the best of our knowledge, the first non-trivial algorithm was published 
in 1971 as part of Dinitz’s celebrated maximum flow algorithm; this algorithm takes 0{mn) 
total update time to solve decremental stSP and stR (see [DinOG] for details). In other words, 
it takes linear time (0(n) time) per update if we delete all m edges; here, n and m are the 
number of nodes and edges, respectively. Independently from Dinitz, Even and Shiloach 
[ES81] developed an algorithm with the same total update time for SSR and SSSP^. This 
algorithm remains the fastest prior to our work. For directed acyclic graphs, Italiano [Ita88] 
gave a decremental algorithm with a total update time of 0{m). For the incremental version 
of the problem, where we only allow insertions of edges, a total update time of 0{m) is 
sufficient in general directed graphs [ItaSG]. For the fully dynamic version of the problem, 
where both insertions and deletions of edges are allowed, Sankowski obtained an algorithm 
with a worst-case running time of per update, resulting in a total update time of 

0(An^'®'^^), where A is the number of updates.^ 

King [Kin99] showed how to extend the algorithm of Even and Shiloach to weighted 
graphs, giving the first decremental single-source shortest path algorithm with total update 

^It was actually published for undirected graphs and it was observed by Henzinger and King [HK95] that 
it can be easily adapted to work for directed graphs. 

^Sankowski’s worst-case update time for the fully dynamic single-source single-sink reachability (stR) 
problem is 
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time 0{mnW), where W is the maximum edge weight (and all edge weights are positive 
integers)^. Using weight rounding [RT87, Coh98, Zwi02, Ber09, MadlO, BerlS] this can 
be turned into a (1 + e)-approximate single-source shortest paths algorithm with total 
update time 0(mnlog lU), where the O(-) notation hides factors polylogarithmic in n. The 
situation is similar for decremental strongly connected components: The fastest decremen- 
tal see algorithms take total update time 0{mn) ([RZ08, Lacl3, Rodl3]). Thus many 
researchers in the field have asked whether the 0{mn) total update time for the decremental 
setting can be improved upon for these problems while keeping the query time constant or 
polylogarithmic [Kin08, Lacl3, RZ08]. 

Our Results. We improve the previous 0(mn)-time algorithms for decremental SSR, 
approximate SSSP, and SCC in directed graphs. We also give algorithms for stR and 
approximate stSP. We summarize our results in the following theorem. In Figure 1 we 
compare the running times of our new algorithms with the previous solution by Even and 
Shiloach for different densities of the graph and Table 1 summarizes our results. 

Theorem 1.1. There exist decremental algorithms for reachability and shortest path problems 
in directed graphs with the following total update times: 

• SSR and SCC: 0(min(m^/®n^/^, 

= . 

• StR: 0(min(m^/‘^n^/^, 

• (1 -|- e)-approximate SSSP: 0(min(m^/®n^/^"'“°^^^, 

• (1 -|- e)-approximate stSP: 

The algorithms are correct with high probability against an oblivious adversary and have 
constant query time. The total update time of the SCC algorithms is only in expectation. Our 
algorithms can maintain (1 -|- e)-approximate shortest paths in graphs with positive integer 
edge weights from 1 to W if W < "■ and e > 1/ log'^ n for any constant c. 


Discussion. Our main result are dynamic algorithms with constant query time and 
sublinear update time for single source reachability and (1 -|- e)-approximate single source 
shortest paths in directed graphs undergoing edge deletions. There is some evidence that it 
is hard to generalize our results in the following ways. 

• (All pairs vs. single source) The naive algorithm for computing all-pairs reachability 
(also called transitive closure) in a directed graph takes time 0{mn) even in the static 
setting. To date no combinatorial algorithm (not relying on fast matrix multiplication) 
is known that gives a polynomial improvement over this running time. Thus, unless 
there is a major breakthrough for static transitive closure, we cannot hope for a 
combinatorial algorithm for decremental all-pairs reachability with a total update 
time of o{mn) and small query time. Further, it was recently shown [HKN“'“15] 
that no algorithm (not even a non-combinatorial one) can solve decremental all-pairs 

®The total update time is actually 0{md), where d is the maximum distance, which could be 0{nW). 
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Comparison of Running Times 



Graph density a 


Figure 1: Running times of our decremental reachability and approximate shortest paths 
algorithms as a function of the density of the initial graph in comparison to the Even-Shiloach 
algorithm with total update time 0{mn). A point {a, 15) in this diagram means that for a 
graph with m = 0 (n“) the algorithm has a running time of 


Problems 

Running Time 

o{mn) when 

Fastest when 

SSR 

0{rnJ l^n^l'^) 

Q(j^3/4j.j5/4+o(1) ^ 

-g 777,3/7^^12/7+0(1)^ 

m = o(n^) 
m = uo{n) 

m = 

n < m < rJ 

rJI'o < 777 < 77 I 3/9 
77 I 3/9 < 777 < 772 

stSP 

(5(j^5/4j^l/2) 
q(j.j,j^2/3j^4/3+o(1) ^ 

m = o(n^) 
m = ijj{n) 

m < 

m > 

(1 + e)-approximate SSSP 

Q(j.j.j^7/6j^2/3+o(1)) 

(9(^3/4^5/4+o(l)j 

m = o(n^) 
m = w(n) 

m < 

m > 

(1 + e)-approximate stSP 

Q(j.j.j^5/4^1/2+o(1)) 
q(j.j.j^2/3j^4/3+o(1) ^ 

m = o(n^) 
m = aj(n) 

m < 

m > 


Table 1: Summary of our algorithms. The running time in the second column is the total 
update time. For (1 + e)-approximate SSSP and (1 + e)-approximate stSP we assume that, for 
some constant c, e > 1/ log'^ n and the largest edge weight is at most The third column 

shows roughly where these algorithms give an improvement over the previous 0 (mn)-time 
algorithms. The fourth column indicates the density of the graph for which the algorithm 
gives the fastest running time. 














reachability faster, based on an assumption called Online Matrix-Vector Multiplication 
(OMv). 

• (Worst-case update time) The total update time of our single-source reachability algo¬ 
rithms is o{mn) which gives an amortized update time of o{n) over a sequence of 0 (m) 
deletions. It recently has been shown by Abboud and Vassilevska Williams [AVW14] 
that a combinatorial algorithm with worst-ease update time and query time of o(n^) 
per deletion implies a faster combinatorial algorithm for Boolean matrix multiplication 
and, as has been shown by Vassilevska Williams and Williams [VWWIO], for other 
problems as well. Furthermore, for the more general problem of maintaining the 
number of reachable nodes from a source under deletions (which our algorithms can 
do) a worst-case update and query time of both o{m) falsifies the strong exponential 
time hypothesis. It might therefore not be possible to deamortize our algorithms. 

• (Approximate vs. exaet) Our SSSP algorithms only provide approximate solutions. It 
has been observed by Roditty and Zwick [RZll] that any exact combinatorial SSSP 
algorithm handling edge deletions that has a total update time of o{mn) and small 
query time implies a faster combinatorial for Boolean matrix multiplication. After 
the preliminary version [HKN14c] of this work appeared, Henzinger et al. [HKN“'“15] 
showed that 0{mn) is essentially the best possible total update time for maintaining 
exact distances under the assumption that there is no “truly subcubic” algorithm for a 
problem called online Boolean matrix-vector multiplication. Both of these hardness 
results apply even for unweighted undirected graphs. Thus, approximation might be 
necessary to break the 0{mn) barrier. 

Organization. We first introduce the notation and basic concepts shared by all our 
algorithms in Section 2. To provide some intuition for our approach we sketch two basic 
s-t reachability algorithms in Section 3. We give the reductions for extending our results 
to SSSP (and SSR) and SCC in Sections 4.1 and 4.2, respectively. In Section 5 we give a 
hierarchical s-t reachability algorithm, which in turn gives algorithms for SSR and SCC. In 
Section 6 we extend this idea to (1 -|- e)-approximate stSP and SSSP. Finally, in Section 7 we 
give an SSR algorithm with improved total update time for dense graphs. We conclude the 
paper by discussing open problems. 

2 Preliminaries 

2.1 Notation 

We are given a directed graph G that might be weighted or unweighted. The graph undergoes 
a sequence of updates, which might be edge deletions or, for weighted graphs, edge weight 
increases. This is called the decremental setting. We denote by A(G) the number of updates 
in G, i.e., the number of edge deletions and edge weight increases. We say that a node v is 
reachable from u (or u can reach v) if there is a path from u to u in G. The distance dcix, y) 
of a node x to a node y in G is the weight of the shortest path, i.e., the minimum-weight 
path, from x to y in G. If there is no path from x to y in G we set dcix, y) = oo. 

Let G = (V, E) be a weighted directed graph, where V is the set of nodes of G and E is 
the set of edges of G. We denote by n the number of nodes of G and by m the number of 
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edges of G before the first update. We denote the weight of an edge {u,v) in G by wg{u,v). 
In the weighted case we consider positive integer edge weights and denote the maximum 
allowed edge weight by W. We assume that W < for some constant c. In unweighted 

graphs we think of every edge weight as equal to 1. Note that for unweighted graphs we 
have A(G) < m and for weighted graphs we have A(G) < mW. 

For every path vr = {vo,V 2 , ■ ■ ■, vfi) we denote its weight in G by w{'k., G) = X]o<i<fc-i 
and its size in terms of number of edges (also called hops) by | 7 r| = fc. We say that the length 
of a path is its weight, but, to avoid ambiguity, we reserve this notion for unweighted graphs 
where the length of a path is equal to its number of edges. For every integer h > 1 and all 
nodes x and y in a directed graph G, the h-hops distance dQ{x,y) is the minimum weight of 
all paths from x to y in G consisting of at most h edges. Note that dciu, v) = dc(u, v). 

For every graph G = {V, E) and every subset of nodes U CV, we define E[U] = E 
and denote by G[U] = {U, E[U]) the subgraph of G induced by U. We set the weight of 
every edge in G[U] equal to its weight in G. For sets of nodes U <ZV and U' C V we define 
E[U, U'] = E n {U X U'), i.e., E[U, U'] is the set of edges (u, v) G E such that u G U and 
V G U'. 

Our algorithms allow problem-specific update and query operations. The algorithms 
we design will have constant query time and we will compare them by their total update 
time over all deletions and, in case of weighted graphs, edge weight increases. They are 
randomized and correct with high probability against an oblivious adversary [BDBK+94] 
who chooses its sequence of updates and queries before the algorithm starts.^ We say that 
an event happens with high probability (whp) if it happens with probability at least 1 — 
for some constant c. We assume that arithmetic operations on integers can be performed 
in constant time. Note that the total update time of decremental algorithms for weighted 
graphs will be il(A(G)), the number of edge deletions and edge weight increases in G. This 
is unavoidable as the algorithm has to read every update in the graph. In unweighted 
graphs the dependence on A(G) does not have to be stated explicitly because there we have 
A(G) < m and the total update time is n(m) for reading the initial graph anyway. We use 
O-notation to hide logarithmic factors, i.e., we write 0{t{m,n,W)) as an abbreviation for 
0{t{m,n,W)log^ n) when c is a constant. Similarly, we use the notation 0{t(m,n)) as an 
abbreviation for 0{t{m, n) Whenever one of our algorithms returns an a-approximate 

result, the approximation factor will be of the form a = 1 -|- e where 0 < e < 1. In this paper 
we assume that e > 1/ log'^ n for some constant c. 

2.2 Problem Descriptions 

Our goal is to design efficient dynamic algorithms for the following problems. 

Definition 2.1 (stR). A decremental single-source single-sink reachability (stR) algorithm 
for a directed graph G undergoing edge deletions, a source node s, and a sink node t, maintains 
the information whether s can reach t. It supports the following operations: 

• Delete(?x, v): Delete the edge {u,v) from G. 

• Query().‘ Return ‘yes’ if s can reach t and ‘no’ otherwise. 

^In particular, the oblivious adversary does not know the random choices of our algorithm an cannot infer 
them by choosing its updates and queries adaptively. 
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Definition 2.2 (SSR). A decremental single-source reachability (SSR) algorithm for a 
directed graph G undergoing edge deletions and a source node s, maintains the set of nodes 
reachable from s in G. It supports the following operations: 

• DELETE(ii, n); Delete the edge {u,v) from G. 

• QuERY(ti); Return ‘yes’ if s can reach v and ‘no’ otherwise. 

Definition 2.3 (SCC). A decremental strongly connected components (SCC) algorithm for 
a directed graph G undergoing edge deletions maintains a set of IDs of the strongly connected 
components of G and, for every node v, the ID of the strongly connected component that 
contains v. It supports the following operations: 

• DELETE(ii, n); Delete the edge {u,v) from G. 

• QuERY(n); Return the ID of the strongly connected component that contains v. 

Definition 2.4 (a-approximate stSP). A decremental a-approximate single-source single¬ 
sink shortest path (stSP) algorithm for a weighted directed graph G undergoing edge deletions 
and edge weight increases, a source node s, and a sink node t, maintains a distance estimate 
5{s,t) such that dcis,!) < 5{s,t) < adG{s,t). It supports the following operations: 

• DELETE(ii, n); Delete the edge {u,v) from G. 

• Increase (u, v, x): Increase the weight of the edge {u,v) to x. 

• Query(); Return the a-approximate distance estimate S{s,t). 

Definition 2.5 (a-approximate SSSP). A decremental a-approximate single-source shortest 
paths (SSSP) algorithm for a weighted directed graph G undergoing edge deletions and edge 
weight increases and a source node s, maintains, for every node v, a distance estimate S(s,v) 
such that dcisjv) < 5{s,v) < adG{s,v). It supports the following operations: 

• DELETE(ii, u); Delete the edge {u,v) from G. 

• Increase (u, v, x): Increase the weight of the edge {u,v) to x. 

• QuERY(n); Return the a-approximate distance estimate S(s,v). 

2.3 Definitions and Basic Properties 

In the following we introduce the basic concepts shared by all our algorithms. The first new 
concept used by our algorithms is the path union of a pair of nodes. 

Definition 2.6 (Path Union). For every directed graph G, every D > 1, and all pairs of 
nodes x and y of G, the path union V{x, y, D,G) CV is the set containing all nodes that lie 
on some path tt from x to y in G of weight wg{t^) < D. 

Note that if the shortest path from x to y in G has weight at most D, then the subgraph 
of G induced by the path union V{x, y, D, G) contains this shortest path. Thus, instead of 
finding this shortest path in G directly we can also find it in this potentially smaller subgraph 
of G. In our algorithms we will be able to bound the sizes of the path union subgraphs, 
which makes them very useful. Observe that for a fixed value of D, the path union can be 
computed in nearly linear time: 



Lemma 2.7 (Linear-Time Path Union). For every directed graph G, every D > 1 and all 
pairs of nodes x and y of G, we have V{x,y,D,G) = {n G U | dG{x,v) + dG{v,y) < D}. 
We can compute this set in time 0{m) in weighted graphs and 0{m) in unweighted graphs, 
respectively. 

Proof. Clearly, if dcix^v) + dGiy,y) < D, then the concatenation of the shortest path from 
X to v and the shortest path from n to y is a path from x to y of weight at most D and thus 
v G V{x, y, D, G). Conversely, if n G V{x, y, D, G), then there is a shortest path vr from x to 
y containing v of weight at most D. Let tti and 7r2 be the subpaths of tt from x to u and 
from v to y, respectively. Then dG{x, v) -|- dG{v, y) < rcG(vri) -|- wg{tt 2 ) = wg{tt) < D. 

Using Dijkstra’s algorithm we compute dG{x, v) and dG{v, y) for every node u G U in 
time 0{m). Afterwards, we iterate over all nodes and check for every node v whether 
dG^x, v) + dG{v, y) < D in total time 0{n). In unweighted graphs we can compute dG{x, v) 
and dG{v, y) for every node x G U in time 0{m) by performing breadth-first search (BPS). □ 

We further observe in the following that path unions are monotone in the decremental 
setting, i.e., no nodes will ever be added to any path union for a fixed value of D while the 
graph undergoes deletions and weight increases This means that once we have computed the 
path union we can also use it for future versions of the graph, as long as we do not want to 
consider larger weights of the paths. 

Lemma 2.8. Let G be a directed graph and let G' he the result of deleting some edge from 
G or increasing the weight of some edge from G. Then for every D > 1 and all pairs of 
nodes x and y of G, V{x, y, h, G') C V{x, y, h, G). 

Proof. Let x G 'P(x, y, h, G'), i.e., there is a path tt from x to y in G' of weight at most D. 
All edges of vr are also contained in G and thus there is a path from x to y in G' of weight 
'Wg{'^) = wg'{t^) < D. Thus, every node of tt is contained in 'P(x, y, h, G) which implies that 
all edges of vr are contained in G[V{x,y,h,G)]. Therefore WG['p(x,y,h,G)]i'^) = wg{tt) < D 
and thus every node on vr, and in particular x, is contained in V{x, y, h,G). □ 

The last property of path unions we will use repeatedly is that we can update them by 
“computing the path union of the path union”. This rebuilding of path unions will be useful 
later to restrict the overlap of path unions of different pairs of nodes. 

Lemma 2.9. For every directed graph G, every D >1, every pair of nodes x and y, and every 
set of nodes Q such that V{x, y, D,G) C Q CV we have V{x, y, D, G[(5]) = V{x, y, D, G). 

Proof. Let x G V{x, y, D, G[Q]), which means that x lies on a path vr from x to y of weight 
at most D. As G[Q] is a subgraph of G, this path is also contained in G (with the same 
weight) and thus x G 'P(x, y, D, G). 

Now let X G 'P(x, y, D, G), which means that x lies on a path tt from x to y of weight 
at most D. By the assumption V{x,y, D,G) C Q every node v' of vr is contained in Q and 
thus TT is contained in G[Q\ (and has the same weight as in G). As x lies on a path from x 
to y of weight at most D in G, we have x G 'P(x, y, D, G[(5]). □ 

We now introduce two concepts that are used frequently, also in other dynamic algorithms. 
The first such concept is the following algorithmic primitive for maintaining shortest paths 
trees up to bounded depth. 
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Theorem 2.10 (Even-Shiloach tree [ES81, HK95, Kin99]). There is a decremental algorithm, 
called Even-Shiloach tree (short: ES-tree), that, given a weighted directed graph G undergoing 
edge deletions and edge weight increases with positive integer edge weights, a source node 
s, and a parameter D > 1, maintains a shortest paths tree from s and the corresponding 
distances up to depth D with total update time 0{mD + A{G)), i.e., the algorithm maintains 
dais, v) and the parent of v in the shortest paths tree for every node v such that dais, v) < D. 
By reversing the edges of G it can also maintain the distance from v to s for every node v in 
the same time. 


Note that if we want to use the ES-tree to maintain a full shortest paths tree (i.e., 
containing all shortest paths from the source), then we have to set D equal to the maximum 
(finite) distance from the source. If the graph is unweighted, then it is sufficient to set 
D = n — 1 and we can thus maintain a full shortest paths tree with total update time 0{mn). 
If the graph is weighted and W is the maximum edge weight, then the maximum distance 
might be as large as (n — 1)W and the algorithm then takes time 0{mnW), which is less 
efficient than recomputation from scratch even if W = n. Similarly, if we want to maintain a 
shortest paths tree up to h < n — 1 hops (containing all shortest paths with at most h edges), 
the algorithm above takes time 0{mh) in unweighted graphs and 0{mhW) in weighted 
graphs. Using a scaling technique [Ber09, MadlO, BerlS], the running time in weighted 
graphs can be reduced if we allow approximation: we can maintain a (I -|- e)-approximate 
shortest paths tree containing all shortest paths up to h hops in total time 0{mn\ogW/e). 

The second concept used repeatedly in our algorithms is sampling nodes or edges at 
random. All our randomized algorithms will be correct whp against an oblivious adversary 
who fixes its sequence of updates and queries before the algorithm is initialized, revealing its 
choices to the algorithm one after the other. It is well-known, and exploited by many other 
algorithms for dynamic shortest paths and reachability, that by sampling a set of nodes with 
a sufficiently large probability we can guarantee that certain sets of nodes contain at least 
one of the sampled nodes. To the best of our knowledge, the first use of this technique in 
graph algorithms goes back to Ullman and Yannakakis [UY91]. 

Lemma 2.11. Let a> 1, let T be a set of size t and let Si, S 2 , ■ ■ ■, Sk be subsets of T of 
size at least q. Let Lf be a subset of T that was obtained by choo,sing each element of T 
independently with probability p = min(x/g, 1) where x = aln{kt) -|- 1. Then, with high 
probability (whp), i.e., probability at least 1 — l/t“, the following two properties hold: 

1. For every 1 <i <k, the set Si contains a node ofU, i.e., SifiLff^f). 


2. \U\ ^“ixt/q = 0{at\ii{kt)/q). 

Proof. For every 1 < z < /c let Ej be the event that SiCiU = 0, i.e., that Si contains no 
node of U. Furthermore, let E be the event that there is a set Si that contains no node of 
Lf. Note that E = Ui<i<fc 

We first bound the probability of the event Ei for 1 < i < A:. The size of Si CiLf is 
determined by a Bernoulli trial with success probability p. The probability that IS'j n t/| = 0 
is therefore given by 


PriEi) = {I - pf^\ < {I - pY 


1 - 


a In (kt) + 1 



1 

gO In (kt)+l 


1 1 

- < - . 

ekH^ ~ 2kH°- 


10 






Here we use the well-known inequality (1 — 1/y)^ < 1/e that holds for every y > 0. Now we 
simply apply the union bound twice and get 

To bound the size of U, let, for every 1 < i < t, Xj be the random variable that is 1 
if the Tth element is sampled and 0 otherwise, i.e., X is 1 with probability p and 0 with 
probability 1 — p. The total size of U is given by the random variable X = 
p be the expectation of X. By a standard Chernoff bound we have 

Vt[X>{1 + 5)p]<^. 

With 5 = 2 we get that the probability that U has size more than 2>p = 3xt/q is at most 

1 _ 1 ^ 1 _ 1 ^1^1 

gAt ^xtjq — gX galnfct-l-1 — ‘2k°‘t^ — 

By a union bound, the probability that one of property 1 (?7 is a hitting set) or property 2 
(?7 is small enough) fails is at most l/t“ and thus the claim follows. □ 

We now sketch one example of how we intend to use Lemma 2.11 for dynamic graphs. 
Consider an unweighted graph G undergoing edge deletions. Suppose that we want the 
following condition to hold for every pair of nodes x and y with probability at least 1 — 1 /n 
in all versions of G (i.e., initially and after each deletion): if dcix^y) > q, then there is 
a shortest path from x to y that contains a node in U. We apply Lemma 2.11 as follows. 
The set T is the set of nodes of G and has size n. The sets Si, S 2 , ■ ■ ■ ,Sk are obtained 
as follows: for every version of G (i.e., after each deletion) and every pair of nodes x and 
y such that dcix, y) > q, we define a set Si that contains the nodes on the first shortest 
path from x to y (for an arbitrary, but fixed order on the paths). As the graph undergoes 
edge deletions, there are at most m < m? versions of the graph. Furthermore, there are 

pairs of nodes. Therefore we have k < such sets. Thus, for the property above 
to hold with probability 1 — 1/n, we simply have to sample each node with probability 
(In kt)/q = i\nn^)/q = (51nn)/y by Lemma 2.11. 

Finally, we explain how to decrease the number of edge weight increases to be considered. 
Remember that in general the number of edge deletions and edge weight increases can only be 
bounded by A(G) = 0{mW). However, by the reduction below, we may assume without loss 
of generality that A(G') = 0(m log W) as we are only interested in (1 -|- e)-approximate 
shortest paths. This for example helps us in reducing the sampling probability for applying 
Lemma 2.11 as the number of versions of G is reduced significantly. 

Lemma 2.12. Given a weighted directed graph G undergoing edge deletions and edge weight 
increases, we can, in constant time per update in G, maintain a weighted directed graph G' 
undergoing edge deletions and edge weight increases such that 

• the number of updates in G' is at most |'log;^^gVF] per edge (which implies X{G') = 
0{mlogi_^_^W)) and 
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• dcix, y) < dc'ix, ?/)<(! + e)dG{x, y) for all nodes x and y. 

Proof. We define G' to be the graph where each edge weight of G is rounded up to the 
nearest power of (1 + e). For every edge {u,v), we have (1 + e)* < w{u,v) < (1 + for 
some i. By setting w\u,v) = [(1 + we get w{u,v) < w'{u,v) < (1 + e)w{u,v). This 

means that the weight of every path tt in G is (1 + e)-approximated by the weight of vr 
in G'. As this is also true for a shortest path, we get dcix^y) < dcfx^y) < (1 + e)dGix,y) 
for all nodes x and y. With every edge weight increase of an edge {u, v) in G we can in 
constant time compute the weight w'{u, v) and check whether it has increased. The number 
of different edge weights in G is [logj^.,.^ W] and therefore the total number of edge deletions 
and edge weight increases in G' is 0(m log^^^ W). Thus, G' is the desired graph. □ 

3 Algorithm Overview for s-t Reachability 

In this section, we illustrate our main ideas by giving simple algorithms for the s-t reachability 
problem (stR). At the heart of all our algorithms is a new way of maintaining reachability or 
distances between some nearby pairs of nodes using small path unions (defined in Section 2.3). 

3.1 Decremental Bounded-Hop Multi-Pair Reachability 

We hrst give an algorithm for solving the following “restricted” reachability problem: We 
are given k pairs of sources and sinks (si,ti), ..., {sk,tk) and a parameter h. We want to 
maintain, for each 1 < i < A;, whether dG{si,ti) < h. 

In this algorithm we use the following set B C F of hubs of size 0(6) (for some 
parameter 6). Let U C F be a set of nodes obtained by sampling each node independently 
with probability ab\n.n/n and let F C i? be a set of edges obtained by sampling each 
edge independently with probability ablnn/m from the initial graph (for a large enough 
constant a). We let B be the set containing U and the endpoints of every edge in F. 
For every pair of nodes u and v we dehne the hub-distance from u to v as dBiu,v) = 
mina;gB dG{u, x) + dG{x, v). 

The hubs are useful in combination with the path unions in the following way. Intu¬ 
itively, instead of directly maintaining whether dG{si,ti) < h, we can maintain whether 
dG[Pisi,ti,h,G)]isi, ti) < h since G[V{si, U, h, G)] (the subgraph of G induced by P(si, U, h, G)) 
contains all Si-U paths of length at most h in G. This could be helpful when 'P{si, U, h, G) 
is much smaller than F. Our hrst key idea is the observation that all {si,ti) pairs for 
which V{si,ti, h, G) is large can use their paths through a small number of hubs to check 
whether dG[P(si,ti,/i,G)]('So Aj) < h (thus the name “hub”). In particular, consider the following 
algorithm. We maintain the distance from each Si to each ti in two ways. 

1. Maintain dsisijti), for each 1 < i < A:, as long as dB{si,ti) < h. 

2. Once dBisi,ti) > h, construct Q{si,ti) = V{si,ti,h,G) and maintain dG[Q(si,ti)]{si,ti) 
up to value h. 

Our algorithm outputs dG{si,ti) < 6, if and only if either dB{si,ti) < 6, or dG[Q(si,ti)]('®o ^ 
h. The correctness is obvious: either dB{si,ti) < h which already implies that dG{si,ti) < h 
or otherwise we maintain dG[Q(si,ti)]isi,ti) which captures all /i-hop Si-U paths. The more 
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important point is the efficiency of maintaining both distances. Our analysis mainly uses 
the following lemma. 

Lemma 3.1 (Either hub-distance or path-union graph is small). With high probability, for 
each i, either dB{si,ti) < h, or G[P{si,ti, h,G)] has at most m.m{m/b,n‘^/b'^) edges. 

Proof Sketch. By the random sampling of the hubs, if G\P{si, ti, h, G)] has more than n/b 
nodes, then one of these nodes, say x, will have been sampled by the algorithm and thus 
contained in B whp (Lemma 2.11). By definition, x lies on some Si-ti path of length at 
most h. Thus, dB{si,ti) < dG{si,x) + dcixjti) < h. Similarly, if G\V{si,ti,h,G)] has more 
than rajb edges, then one of these edges, say (x, y\ will have been sampled by the algorithm 
whp. This means that x will be contained in B, and thus dB{si,ti) < dG{si,x) + dG{x, U) < h. 
As the number of edges of G[P{si,ti, h, G)] can be at most {Visi, U, h, G)p, the stated bound 
follows. □ 

Lemma 3.1 guarantees that when we construct Q{si,ti), its induced subgraph is much 
smaller than G whp, and so it is beneficial to maintain the distance in G[Q{si,ti)] instead of 
G. For each 1 < i < A:, we can maintain dG[Q(si,ti)]('®o by running an ES-tree rooted at 
Si up to distance h in G[Q{si,ti)], which takes time 0(|E[Q(sj,t^)]| • h) (where E[Q{si,ti)] 
denotes the set of edges of G[Q{si,ti)]). As \E[Q{si,ti)] \ < min(m/6, n^/6^), this takes time 
0{khm.m{m/b,'nf/b‘^)). By Lemma 2.7 we can construct each Q{si,ti) in 0{m) time for 
each 1 < i < k, resulting in a total cost of 0{km) for computing the path unions. Finally, 
the time needed for maintaining the hub distances can be analyzed as follows. 

Lemma 3.2 (Maintaining all dB{si,ti)). We can maintain whether dB{si,ti) < h, for all 
1 < i < k, in total time 0{bmh + kbh). 

Proof Sketch. For each hub v G B, we maintain the values of dG{si,v) and dciv, U) up to h 
using ES-trees up to distance h rooted at each hub. This takes 0{hmh) total update time. 
Every time dG{si,v) or dG{v,ti) changes, for some hub v G B and some 1 < z < A;, we 
update the value of dG{si, U) = min„gB(dG('Si, v) + dG{v, ti)), incurring 0(6) time. Since the 
value of each of dG{si,v) and dG{v,ti) can change at most h times, we need Oikbh) time in 
total. □ 

It follows that we can maintain, for all 1 < z < A:, whether dG{si,ti) < h with a total 
update time of 

O^ bmh + kbh -f k^ -f A;/i • min(m/6, rz^/6^) ^ . (1) 

maintain Zi) construct Q{si,ti) maintain (s;, Z^) 

(Lemma 3.2) 

Note that, previously, the fastest way of maintaining, for all 1 < z < A:, whether dG{si,ti) < h 
was to maintain an ES-tree separately for each pair, which takes 0{mhk) time. 

3.2 Decremental s-t Reachability in Dense Graphs 

In the following we are given a source node s and a sink node t and want to maintain 
whether s can reach t. In our algorithm we use a set of centers G C E of size 0(c) (for some 
parameter c) obtained by sampling each node independently with probability aclnn/n (for 
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some large enough constant a). Using these centers, we define the center graph, denoted by 
C, as follows. The nodes of C are the centers in C and for every pair of centers u and v in C, 
there is a directed edge {u, v) in C if and only if dciu, v) < n/c. 

Lemma 3.3 {C preserves s-t reachability). Whp, s can reach t in G if and only if s can 
reach t in C. 

Proof Sketch. Since we can convert any path in C to a path in G, the “if” part is clear. To 
prove the “only if” part, let tt be an s-t path in G. By the random sampling of centers, there 
is a set of centers ci,C 2 , ■ ■ ■ ,Ck on vr such that ci = s, Ck = t, and dcici, Cj+i) < nfc for all i 
whp (Lemma 2.11). The last property implies that the edge (ci,Ci+i) is contained in C for 
all i. Thus s can reach t in C. □ 

Lemma 3.3 implies that to maintain s-t reachability in C it is sufficient to maintain the 
edges of C and to maintain s-t reachability in C. To maintain the edges of C, we simply have 
to maintain the distances between all pairs of centers up to n/c. This can be done using the 
bounded-hop multi-pair reachability algorithm from before with k = (? and h = n/c, where 
we make each pair of centers a source-sink pair. By plugging in these values in Equation (1), 
we obtain a total update time of Oi/bmn/c P ben ePm cn^/P) for this bounded-hop 
multi-pair reachability. To maintain s-t reachability in C, note that C is a dynamic graph 
undergoing only deletions. Thus, we can simply maintain an ES-tree rooted at s in C. As C 
has 0(c) nodes and O(c^) edges this takes time 0{c^). It follows that the total update time 
of this s-t reachability algorithm is 

o(^ hmn/c + bnc + ^rn -\- cr?/t? -|- (2) 

maintain ds(.,.) construct Q(., ■) maintain dcfQj,, (■, ■) maintain dc(s, t) 

(Lemma 3.2) 

By setting^ h = /rnp!'^ and c = (6n)^/^ = , we get a total update time of® 

O This is o{mn) if m = t<;(re^/^). 

3.3 Decremental s-t reachability in Sparse Graphs 

Maintaining s-t reachability in sparse graphs, especially when m = 0(n), needs a slightly 
different approach. Carefully examining the running time of the previous algorithm in 
Equation (2) reveals that we cannot maintain dQ^Q(^.^ .^,y^{u,v) for all pairs of centers u and v 
at all times: this would cost 0(min(cmn/6, cn^/6^))—the third term of Equation (1) —but 
we always need c = u}{b) to keep the first term of Equation (1) to bmn/c = o{mn) and, since 
b < n, cmn/b = uj{mn) and enf /iP = ojirP) = uj{mn). The new strategy is to maintain 
dG[Q{u,v)]{'^-:V) only for some pairs of centers at each time step. 

®Note that we have to make sure that 1 < 6 < n and 1 < c < n. It is easy to check that this is the case 
using the fact that m < n^. 

® Detailed ealeulation: First note that c < < m. So, the term c® is dominated by the term 

mc^. Observe that b = (rd'jraef^'^, thus rdejP = me?. So, the fourth term is the same as the third 
term [me?). Using c = {bnp^^, we have that the first and third terms are the same. Now, the third 
term is mc^ = For the second term, using bn = c®, we have 

bnc = c‘^ = 
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Algorithm. As before, we sample 0(b) hubs 0(c) centers, and maintain dsi', •) between 
all centers which takes Oimhnjc + bnc) time (Lemma 3.2). The algorithm runs in phases. 
At the beginning of each phase z, the algorithm does the following. Compute a BFS-tree 
T on the outgoing edges of every node rooted at the source s in G. If T does not contain 
t, then we know that s cannot reach t anymore and there is nothing to do. Otherwise, 
let L = (ci, C 2 ,..., Cfc) be the list of centers on the shortest path from s to t in T ordered 
increasingly by their distances to s. For simplicity, we let cq = s and Ck+i = t. Note that we 
can assume that 

Vz, dGici,Ci+i) < njc and dG(ci,Ci+ 2 ) > njc. (3) 

The first inequality holds because the centers are obtained from random sampling (Lemma 2.11) 
and the second one holds because, otherwise, we can remove Cj+i from the list L without 
breaking the first inequality. Observe that L would induce an s-t path in the “center graph”. 
Our intention in this phase is to maintain whether s can still reach t using this path; in 
other words, whether dcCi, Cj+i) < nfc for all z. (We start a new phase if this is not the 
case.) We do this using the framework of Section 3.1: For each pair (ci,Cj+i), we know 
that dG{ci,Ci-\-i) < njc when dsici, Ci+i) < njc. Once dB{ci,Ci-\-i) > n/c, we construct 
Q{ci, Cj+i) = V{ci, Ci+x^njc^ G). After each deletion of an edge (zz, v) in this phase, we find 
every index z such that u and v are contained in Q{ci,Ci+i) and, using the static BFS algo¬ 
rithm, check whether (iG[Q(ci,ci+i)](ci, Cj+i) < We start a new phase when dG(ci, Ci+i) 
for some pair of centers (cj,Ci+i) changes from (ZgCcZ) Q+i) ^ to dcici, Ci+i) > njc. 

Running Time Analysis. First, let us bound the number of phases. As for each pair 
(ci,Ci+i) the distance dGCciTi+i) can only become larger than njc at most once, there are 
at most 0{(?) phases. At the beginning of each phase, we have to construct a BFS-tree in 
G, taking 0{m) time and contributing 0{(P‘m) to the total time over all phases. During the 
phase, we have to construct Q(cj, Cj+i) for at most c pairs of center, taking 0{cm) time (by 
Lemma 2.7) and contributing 0{c^m) total time. Moreover, after deleting an edge (u,v), 
we have to update dG[Q(ci,ci+i)](ci, Ci-i-i) every (5(cj,Cj+i) containing both u and u, by 
running a BFS algorithm. This takes 0{\E[Q{ci,Ci+i\\) time for each Q(ci,Ci+i), which is 
0{mlb) whp, by Lemma 3.1. The following lemma implies that every node will be contained 
in only a constant number of such sets (5(ci,Ci+i); so, we need 0(m^/6) time to update 
^^G[Q(ci,ci+i)](ci,Ci+i) over all m deletions. 

Lemma 3.4. For any i and j > i + 3, Q{ci, Ci+i) and Q{cj, Cj+i) are disjoint. 

Proof. First, we claim that dG{ci,Cj+i) > 2n/c. To see this, let G' be the version of the 
graph at the beginning of the current phase. We know that dG'{ci,Ci+ 4 ,) = dG'{ci,Ci+ 2 ) + 
dG'{ci+ 2 ,Ci+ 4 ) > 2 n/c, where the equality holds because at the beginning of the current 
phase (i.e., in G'), every Ci lies on the shortest s-t path, and the inequality holds because of 
Equation (3). Since j -|- 1 > z -|- 4 and both Ci +4 and Cj+i lie on the shortest s-t path in G', 
dc'ici, Cj^i) > dG'(ci) Ci- 1 - 4 ) > 2n/c. The claim follows since the distance between two nodes 
never decreases after edge deletions. 

Now, suppose for the sake of contradiction that there is some node v that is contained 
in both (5(cj,Ci+i) = V{ci,Ci 4 -i,n/c,G) and Q{cj,Cj+i) = 'P(cj, Cj+i, n/c, G). This means 
that V lies in some Cj-Cj+i path and some cj-Cj+i path, each of length at most n/c. This 
means that dG{ci,v) < njc and dG{v.,Cj+i) < njc and therefore we get dG(ci,Cj+i) < 
dG{ci,v) + dG{v,Cj+i) < Injc. This contradicts the lower bound above. □ 
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Thus, the total update time is 


Oy mbnjc + bnc 

maintain dB(-, •) 
(Lemma 3.2) 


+ 



+ 



construct BPS construct Q(ci, c^+i) 
tree in every phase p^ase 


+ 



update da[Q(ci,ci+i)]ici,Ci+i) 


By setting c = and b = /n = ^ we get a running time of O ? 

This is o{mn) if m = o(n^/^). 


4 General Reductions 
4.1 From SSSP to stSP 

In the following we show a reduction from decremental approximate SSSP to decremental 
approximate stSP. The naive way of doing this would be to use n instances of the approximate 
stSP path algorithm, one for every node. We can use much fewer instances by randomly 
sampling the nodes for which we maintain approximatestSP and by using Bernstein’s shortcut 
edges technique [Berl3]. 

Theorem 4.1. Assume we already have the following deeremental algorithm that, given a 
weighted directed graph G undergoing edge deletions and edge weight increases, a source node s, 
and a set of sinks T of size k, maintains, for every sink t ^T, a distance estimate 5{s,f) such 
that dais, t) < d{s, t) < adds, t)) for some a > 1 with constant query time and a total update 
time ofT{k, m, n). Then, for any k < n, there exists a decremental {l+e)a-approximate SSSP 
algorithm with constant query time and total update time 0(T{0{klog (nA)), m, n) + mnA/k, 
where A = |'log;^^g(nIT)], that is correct with high probability against an oblivious adversary. 

Proof. At the initialization we randomly sample each node of G with probability ak In (reA)/re 
(for a large enough constant a). We call the sampled nodes sinks. We use the decremental 
algorithm for a set of sinks from the assumption to maintain a distance estimate S{s,t) for 
every sink t such that dds,t) < 6{s,t) < adds,t). Additionally, we maintain a graph G' 
which consists of the same nodes as G and the following edges: (1) For every edge (re,n) of 
G we add an edge (re, u) of weight WGfu,v) = (1 + e) . (2) For every sink t we 

add an edge {s,t) of weight WG'{s,t) = (1 + (these edges are called shortcut 

edges). We maintain G' by updating the weights of these edges every time there is an edge 
weight increase or deletions in G or the distance estimate d{s,t) of some sink t changes its 
value. On G' we use Bernstein’s decremental (1 + e)-approximate SSSP algorithm [Berl3] 
with source s and hop count h = n/k. This algorithm maintains a distance estimate 5'{s, re) 
for every node re such that dG'{s,v) < 6'{s,v) < (1 + e)dg,(s,re). 

First, observe that dds, re) < dcfs, re) < (1 + e)dds, re) as the new shortcut edges never 
under-estimate the true distances in G. We now claim that dQ,{s,v) < (1 -|- e)adds,v). If 
dds,v) = oo, then the claim is trivially true. Otherwise let vr be the shortest path from s 
to re in G. All edges of this path are also contained in G'. Thus, if tt has at most h edges, 

^ Detailed calculation: First note that the third term [me?) is dominated by the fourth term (mc^). Using 
b = c^n, the first term is the same as the fourth term (mc^). Now, the fourth term is me® = = 

second term, using bn = c^, we have bnc = c® = which is at most 

(using n < m). For the last term, nP jb = nP rP^"^. 
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then dQt{s, v) = dcis, v). If vr has more than h edges, then we know that the set of nodes 
consisting of the last h nodes of tt contains a sink t whp by Lemma 2.11 as the sinks are 
obtained by sampling from the nodes with probability ak In (nA)/n = a In (nA)/ h. Thus, the 
graph G' contains a shortcut edge (s, t) of weight wg'{s, t) < (1 + d)5{s, t) < (1 + e)adG{s, t). 
Now let tt' be the path from s to u that starts with this edge {s,t) and then follows the 
path TT from t to v. Clearly, the path tt' has at most h edges and is contained in G'. As the 
weight of tt' is at most WG'is, f) + (1 + e)dG{t, v) we get 

dG'{s, v) < wg'{s, u) + (1 + e)dG{t, u) < (1 + e)adG{s, t) + dcCL < (1 + e)adG{s, v) . 

Putting everything together, we get that the distance estimate 5'{s,v) fulfills 

dG{s, v) < dG'{s, v) < 5' v) < (1 + e)dGi{s, u) < (1 + e)'^adG{s, u) < (1 + 3e)adGi{s, v) . 

By running the whole algorithm with e' = e/3 we obtain a (1 + e)a-approximation instead of 
a (1 + 3e)Q;-approximation. 

Finally, we argue about the running time. Our running time has two parts. (1) Bernstein’s 
decremental (1 + e)-approximate SSSP algorithm [Berl3] has constant query time and a total 
update time of 0{mhA) (here we also use the fact that WG'{u,t) increases at most A times 
for every sink t). Thus, our decremental algorithm also has constant query time and by our 
choice oi h = n/k, our total update time contains the term mnA/k). (2) Furthermore, we 
have 0{klog (nA)) sinks whp by Lemma 2.11 (and can make the algorithm fail with small 
probability if not). Thus, the total update time of the decremental algorithm for a set of 
sinks from the assumption is T(0(/clog (nA)), m, n). □ 

4.2 From SCC to SSR 

In the following we reduce decremental strongly connected components to decremental single¬ 
source reachability. Our reduction is almost identical to the one of Roditty and Zwick [RZ08] , 
but in order to work in our setting we have to generalize their running time analysis at 
the cost of losing a factor of logn. They show that an 0{mn) algorithm for single-source 
reachability implies an 0{mn) algorithm for strongly connected components. We show that 
in fact o{mn) time for single-source reachability implies o{mn) time for strongly connected 
components. In the following we will often just write “component” instead of “strongly 
connected component”. 

In contrast to the rest of this paper, we will here impose the following technical condition 
on the decremental single-source reachability algorithm: when we update the algorithm 
after the deletion of an edge, the update procedure will return all nodes that were reachable 
before the deletion, but are not reachable anymore after this deletion. Note that all the 
reachability algorithms we present in this paper fulfill this condition. 

Algorithm. The algorithm works as follows. For every component we, uniformly at 
random, choose among its nodes one representative. In an array, we store for every node 
a pointer to the representative of its component. Queries that ask for the component of 
a node v are answered in constant time by returning (the ID of) the representative of u’s 
component. Using the decremental SSR algorithm, we maintain, for every representative w 
of a component G, the sets I{w) and 0{w) containing all nodes that reach w and that can 
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be reached by w, respectively. Note that, for every node v, we have v € C if and only if 
V G I{w) and v G 0{w). After the deletion of an edge (w, v) such that u and v are contained 
in the same component C we check whether C decomposes. This is the case only if, after 
the deletion, u ^ I^w) or v ^ 0{w) (which can be checked with the SSR algorithm of w). 

We now explain the behavior of the algorithm when a component C decomposes into 
the new components Ci,... ,Ck- The algorithm chooses a new random representative Wi 
for every component Ci and starts maintaining the sets I{wj) and 0{wj) using two new 
decremental SSR algorithms. There is one notable exception: If the representative re of C is 
still contained in one of the components Cj, then for this component we do not choose a new 
representative. Instead, Cj reuses w and its SSR algorithms without any re-initialization. 
The key to the efficiency of the algorithm is that a large component Ci has a high probability 
of inheriting the representative from C. 

Before choosing the new representatives we actually have to determine the new com¬ 
ponents Cl, ... ,Ck- We slightly deviate from the original algorithm of Roditty and Zwick 
to make this step more efficient. If re G Cj, then it is not necessary to explicitly compute 
Cj as all nodes in Cj keep their representative w. We only have to explicitly compute 
Cl, , Cj-i, Cj+i,... Cfc. This can be done as follows: Let A denote the set of nodes that 
were contained in I{w) before the deletion of {u, v) and are not contained in I{w) anymore 
after this deletion. Similarly, let B denote the set of nodes that were contained in 0(11)) 
before the deletion and are not contained in 0{w) anymore afterwards. The nodes \\i Au B 
are exactly those nodes of C that are not contained in Cj. Let C denote the subgraph 
of G induced by AU B. Then the components of C are exactly the desired components 
Cl, ..., Cj-i,Cj+i, ... Ck- Note that the sets A and B are returned by the update-procedure 
of the SSR algorithms of w, which allows us to compute AD B. The graph G' can be 
constructed by iterating over all outgoing edges of A U R and the components of C can be 
found using a static SCC algorithm. 

Analysis. The correctness of the algorithm explained above is immediate. For the running 
time we will argue that, up to a log n-factor, it is the same as the running time of the SSR 
algorithm. When the running time of SSR is of the form 0(m“n^), our argument works when 
a > 1 or /3 > 1. To understand the basic idea (for /3 > 1), consider the case that the graph 
decomposes into only two components Ci and C 2 (with ni,n 2 < n and mi, m 2 < m being 
the corresponding number of nodes and edges). We know that one of the two components 
still contains the representative w. For this component we do not have to spawn a new 
decremental SSR algorithm. This is an advantage for large components as they have a high 
probability of containing the representative. The probability of w being contained in Ci is 
ni/n and it is n 2 /n for being contained in C 2 . Thus, the expected cost of the decomposition 
is ©(mfn^ni/n-l-m"n^n 2 /n). We charge this cost to the smaller component, say Ci. As Ci 
has ni nodes, the average cost we charge to every node in Ci is ©(mfnf/n-l- mfnf~^n 2 /n). 
This amounts to an average cost of 0{m^n^~^) per node. As we will charge each node only 
when the size of its component has halved, the total update time is 0(m“n^ log n). 

Theorem 4.2 (From SSR to SCC). If there is a decremental SSR algorithm with constant 
query time and a total update time of 0{nti{m,n) -|- mt 2 {m,n)) such that ti{m,n) > 1, 
t 2 {m,n) > 1, and ti{m,n) and t 2 {m,n) are non-decreasing^ in m and n, then there exists a 

®The technical assumption that ti{m,n) and t 2 {m,n) are non-decreasing in m and n is natural as usually 
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decremental SCC algorithm with constant query time and an expected total update time of 
0{{nti{m, n) + mt 2 {m, n)) log n). 


Proof. Let us first analyze the costs related to the decomposition of a component. Assume 
that the component Cq decomposes into Ci,..., Let n* and mi denote the number of 
nodes and edges in component i, respectively. Let m[ denote the sum of the out-degrees of 
the nodes in Ci in the initial graph (i.e., before the first deletion). Note that is an upper 
bound on m,. Furthermore we have and J2i=i = uiq. 

Assume that the representative w of Cq is contained in Ci after the decomposition. First 
of all, for every j i we have to pay a cost of 0{njti(mj,nj) + mjt 2 {mj,nj)) for initializing 
and updating the decremental SSR algorithm of the new representative of Cj. Second, we 
have to pay for computing the new components. This consists of three steps: (a) computing 
Ad B, (b) computing G', and (c) computing the components of G'. Remember that AD B 
is the union of A, the set of nodes that cannot reach w anymore after deleting {u, v), and B, 
the set of nodes that w cannot reach anymore after deleting {u,v). After deleting {u,v) the 
incoming SSR algorithm of w outputs A and the outgoing SSR algorithm of w outputs B. 
Thus, the cost of computing AdB can be charged to the reachability algorithms of w (which 
have to output A and B anyway). The graph G' is the subgraph of G induced by the nodes in 
AdB. We construct G' by checking, for every node in AUR, which of its outgoing edges stay 
in AUR. This takes time Using Tarjan’s linear time algorithm [Tar72], we can 

compute the strongly connected components of G' in the same running time. As m'- > mj 
the total cost of Cj is 0{njti{mj,nj) + m'jt 2 {mj,nj)) = 0{njti{m,n) + m'jt 2 {m,n)) 

By the random choice of the representatives, the probability that w is contained in Ci is 
UifuQ. Thus, the expected cost of the decomposition of Cq is proportional to 

~ n)) = ^ ^ njti{m, n) — m'jt 2 {m, n) . (4) 


We analyze each of these terms individually. 

Consider first the cost of 0{J2i=ini/nQj2jjLinjti{m,n)). For every pair i,j such that 
i / j we have to pay a cost of 0{ninj{ti{m,n) + n)/no). If n* < nj we charge this cost 
to the component Ci, otherwise we charge it to Cj (i.e., we always charge the cost to the 
smaller component). Note that the component to which we charge the cost has at most no/2 
nodes (otherwise it would not be the smaller one). For a fixed component i, the total charge 
is proportional to 


E 


niUjitiim, n) + n) 
no 


= ni{ti{m,n) + n) 


\ ' ■ / ■ TL ‘ 

———- < nAtiim, n) +n) . 
no 


We share this cost equally among the nodes in Ci and thus charge 0{t{m, n) + n) to every 
node in Ci. Every time we charge a node, the size of its component halves. Thus, every node 
is charged at most log n times and the total update time for the first term in Equation (4) is 
0{nt{m, n) logn). 

Consider now the cost of 0(X)f=i ra))- Note that 

k k k 


E^EK=EKE^=E 


, no ... 

i=\ ^ 


2 = 1 




no 


i=l 


/ no nj 
no 


m. 


the running time of an algorithm does not improve with increasing problem size. 
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We now charge m[{no — ni)/no to every component Ci {1 < i < k). In particular we charge 
(no — ni)jno to every edge (n, v) of the initial graph such that u G Ci. We now argue that in 
this way every edge is charged only O(logn) times, which will imply that the total update 
time for the second term in Equation (4) is 0{mt2{m, n)logn). Consider an edge {u,v) 
and the component containing u. We only charge the edge (n, v) when the component 
containing u decomposes. Let oq denote the initial number of nodes of this component 
and let Op its number of nodes after the p-th decomposition. As argued above, we charge 
(op-i — ap)/ap-i to {u,v) for the p-th decomposition. Thus, for q decompositions we charge 
J2i<p<qi^p-i ~ Op)/®p-i- Now observe that 


E 

i<p<g 


tip_ ^ Qj' 

dp—I 


cip—i dp 1 


dp—\ dp 1 


"< E E 

l<p<g ^=0 


< E E 


“P-I l<p<q i=0 * 


a,p-i ao 1 

= L L (= E 

l<p<q i=dp-\-l i=dq-\-l 


Since oq < n, this harmonic series is bounded by O(logn). 

Finally, we bound the initialization cost. Let Ci,... ,Ck denote the initial components 
and let rii and nii denote the number of nodes and edges of component Ci, respectively. 
The initial components can be computed in time 0{m) with Tarjan’s algorithm [Tar72]. 
Furthermore, each component starts two decremental SSR algorithms and we have to pay 
for the total update time of these algorithms. This time is proportional to 


k k k 

'^{niti{mi, rii) + t 2 {mi,ni)) < ti{m, n)'^ni + t 2 {m, n)'^mi < nti{m, n) + mt 2 {m, n) . 
2 = 1 2=1 2=1 


□ 


5 Single-Source Single-Sink Reachability 

In this section we give an algorithm for maintaining a path from a source node s to a sink 
node t in a directed graph undergoing edge deletions, i.e., we solve the decremental stR 
problem. Using the reduction of Section 4.1, this implies an algorithm for the decremental 
single-source reachability (SSR) problem. 

5.1 Algorithm Description 

Our s-t reachability algorithm has a parameter A: > 1 and for each 1 < z < A: parameters 
bi < n and Ci < n. We determine suitable choices of these parameters in Section 5.3. For 
each 1 < i < A; — 1, our choice will satisfy bi > bi+i and c, > 2ci+i. We also set bk+i = 1, 
Co = n, Cfc+i = 1, and hi = n/ci for all 0 < i < A: -|- 1. Note that this implies /ij+i > 2hi 
for all 1 < z < A;. Intuitively, bi and Cj are roughly the number of z-hubs and z-centers used 
by our algorithm and hi is the hop range of the z-centers. In the algorithm we will often 
consider ordered pairs of centers of the form (x,y). 
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Initialization. At the initialization (i.e., before the hrst deletion), our algorithm determines 
sets of nodes Bi D B 2 ^ ^ Bk and Co 2 Ci 3 • • • 3 Ck+i as follows. For each I < i < k, 

we sample each node of the graph with probability a6j In (A(G))/n and each edge with 
probability aft* In (A(G))/m (for a large enough constant a). The set Bi then consists of the 
sampled nodes and the endpoints of the sampled edges. We set Go = 1^ and G^+i = {s,t}. 
For each 1 < i < fe, we sample each node of the graph with probability acj In (A(G))/n (for 
a large enough constant a). The set Gj then consists of the sampled nodes together with the 
nodes in Gj+i. For every 1 < z < A: we call the nodes in Bi z-hubs and for every 0 < i < k + 1 
we call the nodes in Gj z-centers. Note that, as A(G) < m in unweighted graphs, the number 
of z-hubs is 0(5j) whp and the number of z-centers is 0{ci) by Lemma 2.11. We can make 
the algorithm fail with small probability if these sets are larger. 

Data Structures. Our algorithm uses the following data structures: 

• For every z-hub 2 ; (with 1 < z < fe) an ES-tree up to depth 2hi in G {outgoing tree of 
z) and an ES-tree up to depth 2hi in the reverse graph of G {incoming tree of z) 

• Eor every pair of z-centers {x, y) (with 1 < z < A:) the set of z-hubs from Bi linking x 
to y. 

• For every pair of z-centers {x,y) (with 1 < z < A;) a set of nodes Q{x,y,i), which is 
initially empty. 

• For every pair of z-centers {x, y) (with 0 < z < A:) a list of pairs of (z -|- l)-centers called 
(z -|- l)-parents of {x,y). 

• For every pair of z-centers {x, y) (with l<z<A:-|-l)a list of pairs of (z — l)-centers 
called (z — l)-children of {x,y). 

Maintaining Hub Links. For every 1 < z < A:, we say that an z-hub z links an z-center x 
to an z-center y if dG{x, z) < 2hi and dG{z, y) < 2hi. The hub links of every pair of z-centers 
(x, z/), for all 1 < z < A;, can be maintained as follows. After every edge deletion in the graph, 
we report the deletion to the ES-trees maintained by the hubs. Initially, and after each 
deletion, the level of a node v in the incoming (outgoing) tree of an z-hub 2 ; is at most 2hi if 
and only if dG{v., z) < 2hi {dG{z, v) > 2hi). Thus, we can check whether an z-hub 2 ; links an 
z-center x to an z-center y by summing up the levels of x and y in the incoming and outgoing 
ES-tree of z, respectively. Thus, for every priority z and every pair of z-centers (x, y) we can 
initialize the set of z-hubs linking x to y by iterating over all hubs. We can maintain these 
sets under the edge deletions in G as follows: Every time the level of some z-center x in the 
incoming ES-tree of some z-hub 2 ; exceeds 2/zj, we iterate over all z-centers y and remove 2 : 
from the set of hubs linking x to y. We proceed similarly if the level of some z-center y in 
the outgoing ES-tree of some z hub 2 ; exceeds 2hi. In this way, we can also generate, after 
every edge deletion, a list of pairs of z-centers (x, y) such that x is not linked to y by an 
z-hub anymore (but was linked to an z-hub before the deletion). 

Main Algorithm. Initially, the algorithm computes the shortest path vr from s to t in 
G. It then determines a sequence of A:-centers xi,..., x; on vr. For each 1 < j < ^ — 1, the 
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algorithm now tries to maintain a path from Vj to Wj+i under the edge deletions in G. If 
it fails to do so, it recomputes the shortest path from s to t in G. We call this a refresh 
operation. 

For each 1 < j <1 — 1, the path from Vj to Vj^i is maintained as follows. The first case 
is that there is a A:-hub linking Vj to Uj+i. As long as such a hub exists, we know that there 
is a path from vj to Uj+i. If there is no A:-hub linking Vj to Uj+i anymore, the algorithm 
computes the path union V{vj,Vj^i, 2hk, G). It then recursively maintains the path from Vj 
to Vj+i in G[P{vj,Vj+i,2hk,G)] using {k — l)-centers and {k — l)-hubs. To keep track of this 
recursive hierarchy, we make, for every 1 < j <1 — 1, each pair {vj,Vj^i) an (i — l)-child 
of {x,y) and {x,y) an i-parent of {vj,Vj+i). The end of the recursion is reached in layer 
0, where the set of 0-centers consists of all nodes. The only paths between 0-centers the 
algorithm considers are the edges between them. 

Algorithm 1 shows the pseudocode for these procedures. A crucial ingredient of our 
algorithm is the computation of path unions for pairs of centers. For every 1 < i < k and 
all z-centers {x,y) the path union from the last recomputation is stored in Q{x,y,i). To 
simplify the formulation of the algorithm we set Q{s,t,k + 1) = V (which is consistent with 
our choice of h^+i = n). Thus, the procedure Refresh(s, t, k + 1) simply checks whether 
there exists a path from s to t in G. 

5.2 Correctness 

It is obvious that the algorithm correctly maintains, for all pairs of z-centers (x, y), the set 
of z-hubs linking x to y. Furthermore, the algorithm only stops if it has correctly detected 
that there is no path from s to t in the current graph G. Thus, it remains to show that, as 
long as the algorithm does not stop, there there is a path from s to t. We say that a pair of 
z-centers {x,y) is active if it has at least one (z -f l)-parent. We further define the pair {s,t), 
the only pair of {k + l)-centers, to be active. 

Lemma 5.1. After finishing the delete-operation, for every 0 < i < k 1 and every active 
pair of i-centers {x,y), there is a path from x to y in the current graph G. 

Proof. The proof is by induction on z. If y is linked to x by some z-hub z, we know that in 
G there is a path from x to z as well as a path from z to y. Their concatenation is a path 
from X to y in G. 

Consider now the case that there is no z-hub linking x to y. If z = 0, we know that there 
is an edge from x to y in G as otherwise the pair (x,y) would not be active. If z > 1, let 
{vi,V2), {v 2 ,vfi), ..., {vi-i,vi) with vi = X and vi = y denote the (z -f l)-children of (x,y) 
found at the time of the last refresh of (x, y). Note that all these children are active since 
(x, y) is their z-parent. Thus, by the induction hypothesis, we know that there is a path from 
Vj to Vj^i in G for all 1 < j < I — 1. The concatenation of these paths gives a path from x 
to y in G, as desired. □ 

5.3 Running Time Analysis 

Important Properties. In the running time analysis we will need some properties of the 
path unions computed by our algorithm. 
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Algorithm 1: Algorithm for decremental s-t reachability 

1 Procedure Refresh(x, y, i) 

2 Compute shortest path tt from x to y in G[Q{x, y, i)] 

3 if length of tt > hi then 

4 i = k + 1 then 

5 1^ Stop and output “s cannot reach t” 

6 foreach {i + l)-parent {x',y') of{x,y) do 

7 |_ Refresh(x', y', i + 1) 

8 else 

9 RemoveChildren(x, y, i) 

10 Determine {i — l)-centers xi,..., x; on tt in order of appearance on vr such that 
vi = X, vi = y and, for each I < j < I — 2, Xj+i is the first (i — l)-center 
following Vj on tt at distance at least hi-if2 from Vj. 

II If i = 1 then the 0-centers are all nodes on vr because /iq = 1. 

11 foreach 1 < j < / — 1 do 

12 Make (xj,Xj+i) an (i — l)-child of {x,y) and {x,y) an i-parent of (xj,Xj_|_i) 

13 if Vj+i not linked to Vj by an i-hub then 

14 COMPUTEPATHUNION(xj, Vj+i, i — 1) 

15 ReFRESH(Xj, Xj+1, i — 1) 

16 Procedure ComputePathUnion(x, y, i) 

17 if Q{x,y,i) = 0 then 

18 Let {x',y') be any {i + l)-parent of {x,y) 

19 Q{x, y, i) ^ V{x, y, 2hi, G[Q{x', y', i + 1)]) 

20 else 

21 ^ Q{x, y, i) ^ V{x, y, 2hi, G[Q{x, y, i)]) 

22 Procedure RemoveChildren(x, y, i) 

23 foreach (i — l)-ehild (x',y') of {x,y) do 

24 Remove {x\y') from (i — l)-children of (x,y) and remove (x,y) from i-parents 
of (x',y') 

25 if (x',y') has no i-parents anymore and i > 1 then 

26 RemoveChildren(x', y', i — 1) 

27 Procedure Initialize() 

28 foreach 1 < i < k and all i-eenters (x, y) do Q{x, y,i) ^ 

29 ^('5, t, “t“ 1) ^— R 

30 Refresh (s, t, k+1) 

31 Procedure Delete(u, x) 

32 foreach 1-parent (x,y) of{u,v) do 

33 Refresh(x, y, 1) 

34 For every 1 < i < k and every pair of i-centers (x, y) that is not linked by an i-hub 
anymore: Refresh(x, y, i) 




Lemma 5.2. Let {x',y') be pair of nodes such that dc{x',y') < /ij+i and let {x,y) be a pair 
of nodes such that x and y lie on a shortest path from x to y in G (and x appears before v 
on this shortest path). Then V{x,y,2hi,G) CT’(x',y',2hi^i,G). 

Proof. Let v G V{x,y,2hi,G). We first apply the triangle inequality: 

dcix, v) + dciv, y) < dcix, x) + dcix, v) + dciv, y) + d^y, y) ■ 

Since x and y (in this order) lie on the shortest path from x' to y', we have dcix^x) + 
dciy^y') < dcix',y') and by our assumption we have dcix',y') < /ij+i- Since v G 
7^(x, y, 2hi, G), we have dcix, v) + dciv, y) < 2hi As 2hi < /ij+i by our assumptions on the 
parameters we get 


dcix', v) + dciv, y') < hj+i + 2hi < hi+i + hj+i = 2hi+i 

which implies that v G V{x', y\ 2hi+i, G). □ 

Lemma 5.3. For every 1 <i < k and every pair of i-centers {x,y), if Q{x,y,i) / 0? then 
V{x,y,2hi,G) C Q{x,y,i). 

Proof. The proof is by induction on i. We first argue that it is sufficient to show that after 
the first initialization of Q(x, y, i) in Line 19 we have Q{x, y, i) = V{x, y, 2/ij, G). After the 
initialization, and before Q{x,y,i) is recomputed, Q{x,y,i) might be “outdated” due to 
deletions in G, i.e., not be equal to V{x,y,2hi,G) anymore. However, since V{x,y,2hi,G) 
only “loses” nodes through the deletions in G, it is still the case that V{x,y,2hi,G) C 
Q{x,y,i) C V- Thus, this is also true directly before Q{x,y,i) is recomputed for the first 
time in Line 21. There Q{x,y,i) is updated to T’{x,y,2hi,G[Q{x,y,i)])- By Lemma 2.9 we 
know that V{x, y, 2hi, G[Q{x, y, i)]) = P(x, y, 2/ij, G), i.e., after the recomputation, Q{x., y, i) 
is equal to 'P(x, y, 2hi, G) again. By repeating this argument, we get that V{x, y, 2hi, G) C 
Q{x,y,i) C V. 

We now show that after the first initialization of Q(x, y, i) in Line 19 we have Q(x, y, i) = 
V{x,y,2hi,G). Let {x',y') be the (i + l)-parent of (x,y) used in the initialization. We will 
show that V{x, y, 2hi, G[Q{x'., y',i + l)]) = P(x, y, 2hi, G). For i = k the claim is trivially true 
because (s, f) is the only (A; + l)-parent of (x, y) and Q{s, t, k + 1) = V. For i < k — 1, we know 
by the induction hypothesis that V{x', y', 2/ii+i, G) C Q{x', y', f + 1). We also know that x 
and y lie on the shortest path from x' to y' in G[Q(x', y', i +1)] such that x precedes y on this 
shortest path and the path has length at most hi. As G[Q{x',y'.,i + 1)] is a subgraph of G, 
this path is also contained in G. Thus, d^x', x) + dciy, y') < hi. By Lemma 5.2 this implies 
that V{x,y,2hiG) C 'P(x', y', 2/ij+i, G). It follows that V{x,y,2hi,G) C Q(x',y',i + 1). 
Thus, by Lemma 2.9, we have 'P(x, y, 2/ij, G[Q{x', y', i + 1)]) = V{x, y, 2hi, G) as desired. □ 

Lemma 5.4. Let (x,y) be an active pair of i-centers and let (x',y') be an (i + l)-parent of 
(x, y). Then Q{x, y, i) C Q(x', y', i + 1). 

Proof. Consider the point in time when (x', y') becomes an (z + l)-parent of (x, y). This can 
only happen in the else-branch of Refresh(x, y, i + 1). Then we know that x and y lie on a 
path from x' to y' of length at most /ij+i in G[Q{x', y', i + 1)] such that x precedes y on this 
path and thus dcix', x) + dciy, y') < ^j+i. From Lemma 5.2 it follows that 'P(x, y, 2/ij, G) C 
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V{x',y',2hi^i,G). Furthermore, P(x', y', 2/ij+i, G) C (5(x', y', i + 1) by Lemma 5.3. As the 
path-union of (x,y) is recomputed when (x,y) becomes an i-child of (x',y'), we also have 
Q{x, y, i) = V{x, y, 2hi, G). Therefore we get that Q{x, y, i) C Q{x', y', i + 1). Furthermore, 
every time the algorithm recomputes Q{x', y', i + 1), (x', y') will stop being an (i -|- l)-parent 
of (x, y). It might immediately become in (i -|- l)-parent again and in this case our argument 
above applies again. □ 

Lemma 5.5. For every 1 < i < k + 1 and every pair of i-centers (x,y), we have that if 
dG[Qix,y,i)]{x,y) > hi, then dcix^y) > hi. 

Proof We show that dG{x,y) < hi implies n) ^ If dG{x,y) < hi, then also 

dG['P{x,y, 2 hi,G)]ix, y) < hi as V{x,y,2hi,G) contains all paths from x to y in G of length 
at most 2hi. Since V{x,y,2hi,G) C Q{x,y,i) by Lemma 5.3, we have d(^[Q(a;,j;,i)](x, y) < 
dGir{x,y, 2 hiG)]{x,y)- It follows that dc[Q{x,y,i)]{x,y) < hi. □ 

Lemma 5.6. For every 1 < i < k + 1 and every pair of i-centers (x, y) the graph G[Q{x, y, i)] 
has at most mm{m/bi,n'^/bf) edges whp. 

Proof. As no edges are ever added to G, we only have to argue that the claim is true at the first 
initialization of Q{x, y, i) in Algorithm 1, where Q{x, y, i) is equal to V{x, y, 2hi, G) (see proof 
of Lemma 5.4). We argue that at that time G\P{x,y,2hi,G)] as at most mm.{m/bi,n^/bf) 
edges. Note to this end it is sufficient to show that G['P(x, y, 2hi, G)] has at most n/bi nodes 
and m/bi edges. 

Suppose that V{x,y,2hi,G) contains more than n/bi nodes. Then, by the random 
sampling of i-hubs, one of these nodes, say v, would have been sampled whp while determining 
the set Bi at the initialization by Lemma 2.11, making v an i-hub. If G[P{x,y,2hi,G)] 
contains more than m/bi edges, then one of these edges, say (u, v), would have been sampled 
whp while determining the set Bi at the initialization by Lemma 2.11, making both u and v 
i-hubs. 

In both cases, 'P(x, y, 2/ij, G) contains some i-hub v, for which dG{x, v) -f dG{v, y) < 2hi 
by Lemma 2.7. But this means that x is linked to y by the i-hub v and the algorithm would 
not have executed Algorithm 1, which contradicts our assumption. □ 

Lemma 5.7. Let 1 < i < k and consider a pair of active {i -|- l)-centers {x',y') and their 
i-children {xj,yj)i<j<i (which are active i-centers). Then, for every node v, there are at 
most q = 8 pairs of i-children {xj,yj) of {x',y') such that v G Q{xj,yj,i). 

Proof. We show that at the last time the algorithm has called Refresh(x', y/ i -f 1) (where 
it determined the current i-children of (x^ y')) there are at most q = 8 pairs of i-children 
{xj,yj) of (x',y') such that v G V{xj,yj,2hi,G) as, for each 1 < i < Z, Q{xj,yj,i) (even if 
initialized later) will always be a subset of this set. 

Suppose that v is contained in y > 8 path unions V{xj,yj,2hi,G) of i-children {xj,yj) 
of (x',y'). Let ji,... ,jq he the corresponding indices and assume without loss of generality 
that ji < j 2 ■ ■ ■ < jq- The children of (x', y') all lie on the shortest path from x' to y' and by 
the way we have selected them we even more have dG{xj,yj) >hi/2 for all 1 < j < Z and 
thus 

dG{xji,yj„)> dG{xj,yj)> hi/2 = qhi/2. (5) 
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Furthermore, since v G 'P{xj^,yj^,2hi^G), v lies on a shortest path from Xj-^ to 
of length at most hi in G and thus dcixj^^v) < 2hi. By the same argument we have 
dc{v,yj^) < 2hi. We now have 

dG{xji,yjJ < dG{xj,,v) + dG{v,yj^) < 4hi ( 6 ) 

by the triangle inequality. Observe that the upper bound 6 and the lower bound 5 contradict 
each other for q > 8 and thus q < 8. □ 

Lemma 5.8. For every 1 < i < A: + 1 and every node v, there are at most pairs of 

active i-centers {x,y) such that v G Q{x,y,i). 

Proof The proof is by induction on i. The base case is i = k + 1 and is trivially true 
because the only pair of {k + l)-centers is the pair (s, t). Consider now the case 1 < i < A: 
and fix some node v. Consider a pair (x, y) of active i-centers such that v G Q{x, y, i). Let 
{x',y') be a pair of (i + l)-centers that is an (i -|- l)-parent of {x,y) (such a parent must 
exist because otherwise {x,y) would not be active). Since Q{x,y,i) C Q(x',y',i -f 1) by 
Lemma 5.4, v G Q{x', y', i 4- 1). Thus, v G Q{x, y, i) only if there is an {i -1- l)-parent (x', y') 
such that v G Q{x', y', i -f 1). 

By the induction hypothesis, the number of pairs of active (i-f l)-centers (x', y') such that 
v G Q{x', y', i -f 1) is at most gi^-(*+i)+i = q^-'^. Let (x', y') be such a pair of (i + l)-centers. 
By Lemma 5.7, the number of i-children (x, y) of (x', y') such that v G <3(x, y, i) is at most 
q. Therefore the total number of active pairs of i-centers (x, y) such that v G Q{x, y, i) is at 
most q ■ q^~^ = as desired. □ 

Lemma 5.9. For every 0 < i < k, each active pair of i-centers has at most q^~^ < q^ 
(i + 1)-parents. 

Proof. Consider any active pair of i-centers (x, y) and fix some node v G Q{x,y,i). Let 
the pair of (i + l)-centers (x', y') be an (i + l)-parent of (x, y). By Lemma 5.4 Q{x, y, i) C 
Q{x', y', i -f 1). Thus, v G Q{x', y', i -f 1) for every (i -f l)-parent of (x, y). If (x, y) had more 
than (i -|- l)-parents, then v would be contained in more than path 

unions Q{x',y',i + 1) of (i + l)-parents (x',y'), contradicting Lemma 5.8. □ 

Maintaining Hub Links. For every i-hub, we maintain an incoming and an outgoing 
ES-tree of depth 2hi = 2njci, which takes time 0{mn/ci). As there are 0{hi) hubs of priority 
i, the total time needed for maintaining all these ES-trees is 0{J2i<i<k^i^n/ci). Eor every 
pair of z-centers (x, y), the list of hubs linking x to y is initialized by iterating over all z-hubs. 
As there are 0{ci) z-centers and 0{bi) z-hubs, this takes time 0{J2i<i<kbiCi). Every time 
the level of an z-center x in the ES-tree of an z-hub z exceeds 2/zj, we have to remove z from 
the set of hub links for every possible partner y of x. As this event can occur only once for 
every z-center x and every z-hub z over the course of the algorithm, maintaining the sets of 
linking hubs takes time 0{J2i<i<k 

Computing Path Unions. We now argue about the time needed for maintaining the 
sets Q{x, y, z) for all pairs of z-centers (x, y) and each 1 < z < A. We will show that we 
can pay for this cost by charging 0(min(m/6j+i, to every pair of z-centers and 
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0{q^ min(m/6j+i, n‘^/b‘f_^_l)) to every refresh operation of the form Refresh(x', y', i + 1) for 
some pair of (i + l)-centers {x',y'). Note that when i = k then min(m/6fc+i,n^/6|_(_^) = m 
since we have set bk+i = 1- 

Fixing some pair of i-centers {x,y), we first bound the cost for the first initializa¬ 
tion of Q{x,y,i), as performed in Line 19 of Algorithm 1. There we have to compute 
V{x,y,2hi,G[Q{x',y',i + 1)]), where {x',y') is an (i -|- l)-parent of {x,y). By Lemma 2.7 
this takes time proportional to the number of edges in G[Q{x',y\i -|- 1)]. If i = k, then 
Q{x',y',i + = V (and actually x' = s and y' = t) and thus computing V{x,y,2hi,G) for 
all pairs of fe-centers (x, y) takes time 0{c\m). If 1 < i < fe — 1, then G[Q{x', y', i -|- 1)] has 
0{m.m{m/bi+i,'n?/bf_^_l)) edges by Lemma 5.6. Thus, the first computation of Q{x,y,i) for 
all pairs of Tcenters {x,y) takes time 0{cf m.m{m/bi+i,n'^/b^_^_l)). 

Now consider the cost of computing Q{x, y, i) after it has already been initialized 
for the first time. Let Q'{x,y,i) = V{x,y,2hi,G[Q{x,y,i)]) denote the updated path 
union as computed in Line 21 of Algorithm 1. The cost of computing Q'{x,y,i) is pro¬ 
portional to \E[Q{x,y,i)]\, the number of edges in G[cQ{x,y,i)] before the recomputa¬ 
tion. Since Q'{x,y,i) C Q{x,y,i), we have E[Q{x,y,i)] = E[Q'{x,y,i)] U E[Q{x,y,i)] \ 
E[Q'{x,y,i)]. Note that Q'{x,y,i) is equal to V{x,y,2hi,G). We pay for this cost by 
charging 0{\E[Q{x,y,i)] \ E[Q'(x, y, i)]\) to the pair {x,y) and \E[V{x,y,2hi,G)]\ to the 
refresh operation on the (i -|- l)-parent of (x, y) which causes the recomputation. 

As the initial size of G[Q{x, y, i)] is 0(min(m/6i+i, n^/6?_,_^)) and we only charge edges to 
the pair (x, y) that will never be contained in G[Q{x, y, i)] anymore, the total time charged to 
(x,y) is 0{m.m{m/bi+i,'n?/b‘f_^_l)), which, over all path union computations, results in a cost 
of 0(X]i<i<fc-i cf nim{m/bi+i,n‘^/bf^i) -|- c\m). It remains to bound the total cost charged 
to each refresh operation Refresh(x', y', i-|-l) for some pair of (i-|-l)-centers {x\ y'). Below 
we will then separately analyze the total cost of the refresh operations. During the refresh 
operation we recompute Q{x, y, i) for every z-child (x, y) of (x', y'), for which we have to pay 
0{\E[V{x,y,2hi,G)]\) per child. By Lemma 5.4 Q{x,y,i) C Q{x',y',i + 1) for every i-child 
(x, y) of {x', y'). By Lemma 5.8, each node is contained in at most path unions of i-children 
of {x\y'). Therefore every edge of G[Q{x',y\i + 1)] is contained in the subgraphs of G 
induced by the z-path union of at most q^ z-children of (x', y'). Thus, the total cost charged 
to the refresh operation is 0{q^\E[Q{x',y'+ 1)]|), which is 0((?^ min(m/6i+i, n^/6?_,_^)). 

Cost of Refresh. Excluding the recursive calls, each refresh operation of the form Re- 
fresh(x, y, z), where 1 < i < k + 1 and (x,y) is a pair of z-centers, is dominated by two 
costs: (1) the time needed for computing the shortest path from x to y in G[Q{x,y,i)] and 
(2) the time needed for computing Q{x',y',i — 1) for every (z — l)-child {x',y') of (x,y) if 
z > 2. As G[Q{x,y,i)] has at most min(m/6i+i, n^/6?_)_^) edges by Lemma 5.6, computing 
the shortest path takes time 0{mm(m/bi+i,'n?/b‘f_^_l)). We have argued above that to pay 
for step (2) the time we charge to that particular refresh is 0{q^ mm{m/bi+i,n^/b‘f_^_^)). It 
remains to analyze how often the refresh operation is called. 

We say that a pair of z-centers (x, y) (for 0 < z < fc) causes a refresh if the algorithm 
calls Refresh (x', y', z -I- 1) (in line 7), where {x',y') is an (z -|- l)-parent of (x,y), after 
detecting that the distance from x to y in G[Q(x,y,z)] is more than hi. If a fixed pair of 
z-centers (x, y) causes a refresh after some deletion in the graph, it will do so for each of 
its (z -|- I)-parents. By Lemma 5.9 the number of (z -|- l)-parents is at most q^. Note that 
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each pair of i-centers {x, y) will only cause these refreshes of its parents once. At the 
time the algorithm makes (x,y) the child of some {i + l)-center we have dcix^y) < hi whp 
as by the initial random sampling of i-centers, every shortest path consisting of hi/2 — 1 
edges contains an i-center whp (Lemma 2.11). Furthermore, the pair (x,y) will never cause 
a refresh anymore in the future as the refresh implies that dcix, y) > hi by Lemma 5.5 and 
thus (x, y) will never be active anymore. 

Now whenever we refresh a pair of [i + l)-centers {x',y'), we charge the running time 
of 0{q^ m.m{m/bi+i,n‘^/b‘f_^_i)) to the i-child {x,y) causing the refresh. By the argument 
above, each pair of i-centers will be charged at most q^ times. Thus, the total time 
needed for all refresh operations on pairs of (i + l)-centers over the course of the algorithm 
is 0{q‘^^cf mm{m/bi+i,n^/b‘f_^_^)) if i > 1. For i = 0, we can bound this by 0{q^^m • 
min(m/6i, )) because every node is a 0-center and thus a pair of 0-centers (x,y) can 
only be active if the graph contains the edge {x,y). 

Total Running Time. Putting everything together, the total running time of our algo¬ 
rithm using k layers is 
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We hrst balance the terms to obtain a running time of We achieve this by 

setting the parameters to k = [log log m] and, for every 1 < i < /c. 
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With this choice of the parameters we get (for all 1 < i < /c), bi > 6i+i, c* > 2cj+i, 
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where the last inequality holds due to m < and i > 0. Thus, the total update time is 

/ 5'2*=-3 2(*^+3)-2 \ 
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Now observe that 


^.2^_3 _2 53 1 111 53 1 

j ^ 2 ('=+ 2)_3 j ^ 2 (''+ 2)-3 = 2(^+^'l-3fl^ 2 2(k+2) ^ 2 (*+ 2 )- 37^5 


By our choice of /c = [log log m] we have 

3. 1 _ J_ 1 1 

Ul '^ 2(*+3)—3 <777,2'“ < 777 2'°g'°S"^ = 777log"i = 2 . 

Thus, the running time of our algorithm is With k = [log log m] and 

q = 8 this is Oirr^l'^n^l'^). 

We now balance the terms to obtain a running time of We achieve 

this by setting the parameters io k = [■\/logn/log q\ and, for every 1 < i < /c, 

l)/(3/ii+l) 

^ ^(4i-A:-3)/(6A:+2) 
a = 2*'-*n2V(3fc+l)j^(3fc+l-4i)/(6fc+2) _ 


With this choice of the parameters we get (for all 1 < i < /c), 6* > 6i+i, Ci > 2cj+i, 

= j^4fc/(3fc+l)^(2fc+2)/(3fc+l) 

Cj “ 6i 

— <clm< 22fc774Ai/(3fc+l)^(2fc+2)/(3fc+l) ^ 
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and 


^.^2 < 22fc77(fc+6*-l)/(3fc+l)j„(7fc+5-12i)/(6fc+2) 

= 22fc77(fc-b/(3fc+l)j^(3fc+l)/(6fc+2)^(2fc+2)/(3fc+l)^6j/(3fc+l)^-6i/(3fc+l) 
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where the last inequality holds due to n < m and m < n?. Thus, the total update time is 

6 (k2‘^k^2k^4k/{3k+l)^{2k+2)/{3k+l)'^ ^ 

Now observe that 

2fc+2 4fc 2,4 1 4 4 1 2,4 1 4 2,1. 4 

777 3fc + l 77 3fc+l = 7773 ' 3'3fc+l 773 3'3fe+l < 7773“'"3'3fe+l 773 < 777 3 + fc773 

By our choice of /c = [-yiog n/log gj we have 2^ < q^ < Since k < logn and q = 8 we 

thus obtain a total update time of 

Theorem 5.10. There is a decremental stR algorithm with constant query time and total 

' ii'nnnfp fiTUP 

0(min(m5/"ni/2^m2/377^/3+«(i))) = 0(mn6/7+°(i)) 
that is correct with high probability against an oblivious adversary. 
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5.4 Extension to Single-Sonrce Reachability 

The algorithm above maintains reachability from a single source s to a single sink t. We can 
easily modify the algorithm to maintain reachability for a set of source-sink pairs (s*, h)i<i<p. 
We simply run p instances of the algorithm, each with a different source-sink pair. Note 
however that the algorithm can use the same set of hubs and centers for all instances. Thus, 
the cost of maintaining the ES-trees does not have to be multiplied by p. We therefore get a 
total running time of 
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and 


^.^2 ^ 22'=p{*:+l-3d/(3fc+l)^(fc+6i-l)/(3A:+l)^(7fc+5-12i)/(6fc+2) 

(3fc+l) j^4fe/ (3fc+l)^(2A:+2)/(3fc+l) 

By the same choices of k as in the s-t reachability algorithm above we get a running times 
of and respectively, for maintaining reachability 

between p source-sink pairs. 

Corollary 5.11. There is a decremental algorithm for maintaining reachability of p source- 
sink pairs with constant query time and total update time 

that is correct with high probability against an oblivious adversary. 

Using the reduction of Theorem 4.1 this immediately implies single-source reachability 
algorithm with a total update time of 0 {mf (we balance the terms and 

mn/p by setting p = and (balance and mn/p by 

setting p = 

Corollary 5.12. There is a decremental SSR algorithm with constant query time and total 

iinflnfp fiTinp 

that is correct with high probability against an oblivious adversary. 

Furthermore, the reduction of Theorem 4.2 gives a decremental algorithm for maintaining 
strongly connected components. 

Corollary 5.13. There is a decremental SCC algorithm with constant query time and 
expected total update time 

0(min(m^/6^2/3^^3/4^5/4+o(i))) ^ 

that is correct with high probability against an oblivious adversary. 

6 Approximate Shortest Path 

In the following we assume that G is a weighted graph with integer edge weights from 1 to 
W undergoing edge deletions. We are given some parameter 0 < e < 1 and our goal is to 
maintain (1 -|- e)-approximate shortest paths. Similar to the reachability algorithm, we first 
give an algorithm for maintaining a (1 -|- e)-approximate shortest path from a given source s 
to a given sink t. We then extend it to an algorithm for maintaining (1 -|- e)-approximate 
shortest paths from a fixed source node to all other nodes. We give an algorithm that 
maintains a (1 -|- e) ^^"’'^-approximate shortest path from s to t. Note that by running the 
whole algorithm with e' = cjifTk -|- 2) instead of e we obtain a (1 -|- e)-approximate shortest 
path. Throughout this section we assume that e > 1/log'^n and W < for some 

constant c 
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6.1 Preliminaries 


The first new aspect for approximate shortest paths, compared to reachability, is that we 
maintain, for each 1 < i < fe + 1 and every pair of i-centers (x, y), an index r{x, y, i) such 
that dQ{x,y) > (1 + and, whenever {x,y) is active, dcix^y) < (1 + Thus, 

r(x, y, i) indicates the current range of the distance from x to y. Roughly speaking, the 
algorithm will maintain a sequence of centers and for each pair of consecutive centers (x, y) 
a path from x to y of weight corresponding to the range of the distance from x to y. This 
will be done in a way such that whenever the algorithm cannot find such a path from x to y, 
then the range has increased. By charging this particular increase in the distance from x to 
y, it can afford updating the sequence of centers. As in the case of reachability such paths 
from X to y will either be found via a hub or in a path union graph. 

The second new aspect is that now we maintain approximate shortest paths up to hi hops 
for certain pairs of z-centers, whereas in the reachability algorithm hi was an unweighted 
distance. This motivates the following modihcation of the path union that limits allowed 
paths to a fixed number of hops. 


Definition 6.1. For all nodes x and y of a graph G and all integers h > 1 and D > 1, the 
/i-hop path union V^{x,y,D,G) is the set containing every node that lies on some path vr 
from X to y in G that has |7r| < h edges (hops) and weight w^tt, G) < D. 

Observe that V{x,y, D,G) = V^{x,y, R,G) by Definition 2.6. Dealing with exact 
bounded-hop paths is usually computationally more expensive than dealing with unbounded 
paths, even in the static setting.® The bounded hop path unions will henceforth only be 
used in the analysis of the algorithm. In the algorithm itself we will compute unbounded 
(i.e., n-hop) path unions in graphs with modified edge weights. For every h > 1 and every 
r > 0 we define a graph that has the same nodes and edges as G and in which we round 
the weight of every edge {u,v) to the next multiple of e(l + /h by setting 


WQh,r{u,v) 


WG{u,v)-h'\ e{l + eY 
e(l + e)®- h 


Note that the definition of implicitly depends on e as well. To gain some intuition for 
these weight modifications, observe that for every path vr in G consisting of h' edges we have 


WQh,r{TT) < wg{t^) + h‘ 


/ e(l + er 


h 


Thus, for every path tt of weight wg{t^) > (1 + consisting of at most h edges we have 


< (1 + e)^ + e(l + eY = (1 + eY ^^. 


As additionally wg{u, v) < WQh,r{u, v) for every edge (u, v) we obtain the following guarantees. 


Lemma 6.2. Let h > 1 and r > 0. For all pairs of nodes x and y we have dG{x,y) < 
dQh,r{x, y) and if d^A, y) > (1 + A, then dQh,r{x, y) < {1 + e)''+^. 

®For example, one can compnte the shortest paths among all paths with at most h edges from a sonrce 
node in time 0 {mh) by running the first h iterations of the Bellman-Ford algorithm. This is not efficient 
enough for our purposes, as we would only like to spend nearly linear time. 
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Now observe that /i-hop path unions can be computed “approximately” by computing 
the (unbounded) path union in which can be done in nearly linear time (Lemma 2.7). 
In our case the unbounded path union in will contain the h-hop path union and we 
will later argue that all unbounded path unions computed by our algorithm have small size. 
In our analysis the following lemma will have the same purpose as Lemma 2.9 for the s-t 
reachability algorithm. 

Lemma 6.3. Let {x,y) be a pair of nodes, let h >1, r > 0, and a > 1, and let Q V 
be a set of nodes such that V^{x,y,a{l + e)”,G) C Q. Then V^{x,y,a{l + e)”,G) C 
V{x,y,a{l + eY+\G^^^[Q]). 

Proof. Let v £ V^{x, y, a(l + e)”, G), i.e., there is a path vr from x to y in G with at most h 
edges and weight at most a(l + e)” containing v. The weight of vr in G^’” is 

WQh,T (vr) < wg{u, v) + h - ^ < a(l + e)” + e(l + e)” < a(l + e)”+^. 

By our assumption all nodes of vr are contained in Q and thus vr is a path in &’^[Q] of 
weight at most q;( 1 + e)”"''^. Thus, all nodes of vr, including v, are contained in T’{x, y, a(l + 
e)”+i,G^’”[Q]). □ 

The advantage of the graph G^’” is that its edge weights are multiples of e(l + ef jh. 
Thus, after scaling down the edge weights by a factor of h/(e(l + e)”) we still have integer 
weights. This observation can be used to speed up the pseudopolynomial algorithm of Even 
and Shiloach at the cost of a (1 + e)-approximation [Ber09, MadlO, Berl3]. 

6.2 Algorithm Description 

Our approximate s-t shortest path algorithm has a parameter k > 1 and for each 1 < i < k 
parameters bi < n and ct < n. We also set bk+i = 1, cq = n, Ck+i = 1, and hi = nfci 
for all 0 < i < A; + 1. Our algorithm determines sets of nodes Bi L B 2 ^ Bk and 

Co D Cl D ■ ■ ■ D Gfc_|_i by random sampling as described in Section 5.1 for the s-t reachability 
algorithm, i.e., i-hubs are obtained by sampling nodes with probability a 6 i In (A(G))/n 
and edges with probability aft* In (A(G))/m and i-centers are obtained by sampling nodes 
with probability acj In (A(G))/n. Note that, as we may assume A(G) < m\ogij^,,W by 
Lemma 2.12 and we further assume that e > 1/log'^n and W < 2^°® for some constant c, 
the number of i-hubs is 0{bi) whp and the number of i-centers is 0{ci) whp by Lemma 2.11. 
We can make the algorithm fail with small probability if these sets are larger. For each 
1 < i <k we order the set of i-hubs in an arbitrary way such that Bi = { 61 , 62 , ■ ■ ■, 

Our algorithm uses the following variables and data structures: 

• An estimate 6{s, t) of the distance from s to t in G 

• For every pair of z-centers (x, y) (with 1 < i < A:) an index r(x, y, i) such that 
0 < r(x, y, i) < logi_|_g(nlT) for the current range of the distance from x to y. 

• For every i-hub z (with 1 < i < A:) and every 0 < r < [logi_,_£(nlT)J an ES-tree 

up to depth (1 + e)^+*+i in CY^^Bx (^outgoing tree of z) and an ES-tree up to depth 
(1 -f g)'’+®+i in the reverse graph of (^incoming tree of z). 
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• For every pair of i-centers {x,y) (with 1 < i < k) an index l(x,y,i) such that 
1 < l{x,y,i) < \Bi\ + 1 that either gives the index l{x,y,i) of the i-hub bn^x,y,i) ^ 
that links x to y for the current value of r{x, y, i) or, if no such i-hub exists, is set to 
l{x,y,i) = \Bi\ + 1. 

• For every pair of i-centers {x,y) (with 1 < i < k) a set of nodes Q{x,y,i), which is 
initially empty. 

• For every pair of i-centers (x, y) (with 0 < i < k) a list of pairs of (i + l)-centers called 
{i + l)-parents of (x, y). 

• For every pair of i-centers (x, y) (with l<i<A: + l)a list of pairs of (i — l)-centers 
called {i — l)-children of (x,y). 

• For every 1 < i < k, a list Ai of z-centers called active i-centers. 


Besides the generalizations for weighted graphs introduced in Section 6.1, our approximate 
s-t shortest path algorithm also deviates from the s-t reachability algorithm in the way it 
maintains the hub links. The main difference to before is that now we maintain the hub 
links only between pairs of active centers. In the new algorithm we say, for 1 < z < /c, that 
an z-hub z € Bi links an z-center x to an z-center y if 


dQShi/e,r(x,y,i){x, Z) d^Sfn / e,r(x ,y ,i) {z, y) ^ (l + c) 


r{x,y,i)+ 2 i+l 


The algorithm can check whether z links x to y by looking up dQS,hje,r{x,y,i){x, z) and 
dQShje,r(x,y,i) {z, y) iu the incoming and outgoing ES-tree of z, respectively. For every 1 < z < /c, 
and every pair of active z-centers x and y, the algorithm uses the index Z(x, y, z) to maintain 
an z-hub bi(^x,y,i) that links x to y. If the current z-hub bi(^y.,y,i) does not link x to y anymore, 
the algorithm increases l{x,y,i) until either such an z-hub is found, or l(x,y,i) = \Bi\ + 1. 
As soon as l{x,y,i) = \Bi\ -|- 1, the algorithm computes the path union between x and y 
and maintains an approximate shortest path from x to y in the subgraph induced by the 
path union. If the algorithm does not hnd such a path anymore it increases the value of 
r(x,y,z) and sets Z(x,y,z) to 1 again, i.e., making the hrst z-hub the candidate for linking 
X to y for the new value of r(x,y,z). The generalization of the s-t reachability algorithm 
is now straightforward and the pseudocode can be found in Algorithm 2. Before the hrst 
update the algorithms executes Procedure Initialize. After each increase of the weight of 
an edge {u,v) or after its deletion the algorithm calls UPDATE(zt, v). 


6.3 Correctness 

To establish the correctness of the algorithm we will show that 

dcis, t) < 6{s, t) < (1 -h ef^^^dcis, t). 

By running the whole algorithm with e' = cjiAk -|- 2) (instead of e), we then obtain a 
(1 -|- e)-approximation (see Lemma 6.7 below). 

We will achieve this by showing that with the value of r(x, y, z) of a pair of z-centers (x, y) 
the algorithm keeps track of the range of the distance from x to y. In particular, we will show 
that dg (x, y) > (l-|-e)^(*’2^’d for every pair of z-centers (x, y) and doix, y) < (l + e)’’(*’ 2 ^’*)+ 2 *+i 
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Algorithm 2: Algorithm for decremental approximate s-t shortest path 

1 Procedure Refresh(x, y, i) 

RemoveChildren(x, y, i) 
if i = k + 1 then 

I Compute shortest path vr from s to f in G and set 5{s,t) (1 + e)‘^^^^dG{s,t) 

else 

Compute shortest path vr from re to y in &*’'^^^’y’^^[Q{x,y,i){x,y)] 
if y)] > (1 + then // Range has increased 

r{x,y,i) ^ r{x,y,i) + 1 
l{x,y,i) ^ 1 

foreach {i + l)-parent {x',y') of (x,y) do Refresh(x', y' , i + 1) 

break 

Determine (i — l)-centers vi,... ,vi in order of appearance on vr 
foreach 1 < j < Z — 1 do 
Add {vj,Vj+i) to Aj_i 

Make {vj,Vj+i) an {i — l)-child of (x,y) and (x,y) an Z-parent of {vj,Vj-^-l) 

if i > 1 then 

UPDATEHuBLiNKS(nj, Vj+i, 1 - 1, r{vj,Vj+i,i - 1)) 
if li{vj,Vj^i) = \Bi\ + 1 then 
Q{vj,Vj+i,i- 1) ^ 
r{vj,Vj+l, (1 + 

REFRESH(nj, Vj+i, i - 1) 

21 Procedure UpdateHubLinks(x, y, i, r) 

22 while l(x,y,i) < \Bi\ and 

(^G^hi/e,r(x,y,i){x,bpx,y,i)) dQShi/e,r(x,y,i){byx,y,i)-iy) ^ ^ dO 

23 l{x,y,i) ^ l{x,y,i)+ 1 

24 if l{x,y,i) = \Bi\ + 1 then 

25 Let (x',y') be any (z + l)-parent of (x,y) 

26 [ Q{x,y,i) + + 

27 Procedure RemoveChildren(x, y, i) 

28 foreach (i — l)-child (x',y') of {x,y) do 

29 Remove (x', y') from (i — l)-children of (x, y) and (x, y) from z-parents of (x', y') 

30 if (x',y') has no i-parents anymore then 

31 Remove {x',y') from Aj_i 

32 if z > 1 then RemoveChildren(x', y' , z — 1) 

33 Procedure Initialize() 

34 foreach 1 < z < A: and all i-centers (x, y) do r{x,y,i) •(— 0, l{x,y,i) <— 1, 
Q{x,y,i) ^ 0 

35 foreach l<z</cdo Ai 

36 Refresh (s, t, k + 1) 

37 Procedure Update (zz, v) // Deletion or weight increase involving {u,v) 

38 foreach 1-parent (x,y) of{u,v) do Refresh(x, y, 1) 

39 foreach 1 < z < A: and (x, y) G Ai do 

40 I UpdateHubLinks(x, y, i, r{x,y,i)) 
if l{x,y,i) = \Bi\ + 1 then Refresh(x, y, z) 
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for every pair of active i-centers (x, y) G Ai. Both inequalities are obtained by a “bottom up” 
analysis of the algorithm. In order to prove the lower bound on y) we have to show 

that the subgraph Q{x, y, i) computed by the algorithm indeed contains the corresponding 
path union. In fact those two invariants rely on each other and we will prove them mutually. 

Lemma 6.4. Algorithm 2 correctly maintains the following invariants whp: 

(11) For every 1 <i < k and every pair of i-centers (x,y), 

4“(x,y)>(l + er("’^’*). 

(12) For every 1 <i < k and every pair of active i-centers (x,y) G Ai, 

y, (1 + c y, i ). 


Proof Both invariants trivially hold directly before the first refresh in Algorithm 2 of the 
initialization is called because we set r(x, y,i) = 0 for all 1 < i < A: and no pairs of centers 
are active yet. 

We first give a proof of Invariant (II). As distances are non-decreasing in G we only 
have to argue that the invariant still holds after each time the algorithm changes the value 
of r(x, y, i). The only place where r(x, y, i) is changed is in Line 8 (where it is increased by 
1) and if this line is reached the condition 


j)](x, y) > (1 -I- e) 


r{x,y,i)-{-2 


holds. 

Suppose that in G there is a path vr from x to y with at most hi edges and weight less 
than (1 -|- xhe weight of vr in is 


WQh^,r(^,y,i){7r) < wg{tt) -h hi - - - 

hi 

< (1 + Y^{x,y,i)+l Y^ix,y,i) 

< (1 + ^Y{x,y,i)+2 

< (1 _|_ _ 


By Invariant (12) all nodes of vr are contained in Q{x,y,i). Therefore 

dQln,r{x,y,i)yQ(^^^y^iy^{x,y) < WQhi,r(x,y,i){lT) < (iH-c) ^ 

which contradicts the assumption that the algorithm has reached Line 8. Thus, the path vr 
does not exist and it follows that 4(x,y) > (1 + as desired. Therefore Invariant 

(II) still holds when the algorithm increases the value of r(x,y, i) by 1. 

We now give a proof of Invariant (12). The proof is by induction on i. There are two 
places in the algorithm where the value of Q(x, y, i) is changed: Line 19 and Line 26. Note 
that r{x,y,i) is non-decreasing and l{x,y,i) decreases only if r{x,y,i) increases. Therefore, 
for a fixed value of r(x, y, i), the algorithm will “recompute” Q(x, y, i) in Line 19 only if it 
has “initialized” it in Line 26 before. 
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Consider first the case that Q{x,y,i) is “initialized” (Line 26) and let {x',y') be the 
{i + l)-parent of (x, y) used in this computation. If i < k — 1, then by the induction hypothesis 
we know that 


pSh,+i/e^x', y\ (1 + e)”(*'-2/'-*+i)+2(*+i)^ c Q{x\ y\ i + l). 


(7) 


li i = k, then this inclusion also holds because in that case x' = s and y' = t (as s and t 
are the only [k + l)-centers) and we have set Q{s, t,k + 1) = V. We will show below that, 
at the last time the algorithm has called Refresh(x', y', i + l) (where it sets (x,y) as an 
i-child of {x',y') and {x\y') as an (i + l)-parent of (x,y)), we have 

iP®'*‘/^(x, y, (1 + G) C ip8^'+iA(x', y', (1 + ,i+i)+ 2 {i+i) 


and thus (together with (7)) 


y, (1 + e)dw)+2i^ ^ ^ _ 


( 8 ) 


As distances in G are non-decreasing, (8) will still hold at the time the algorithm sets 


Q(x, y, i) = R(x, y, (1 + ^ 


(9) 


according to Line 26. By Lemma 6.3, (8) will imply that 


iP®^^/"(x,y,(l + e)”(^’^’*)+2^G) 

C R(x,y, (1 + e)”("’^’*)+ 2 i+i^g 8 +A,Kx,,,)[g(^/^y/^ . ^ ^ g(x,y,i) 


as desired. Note again that, since distances in G are nondecreasing, the inclusion R®^*/'^(x, y, (1+ 
G) C Q(x, y, i) will then continue to hold until Q(x, y, i) is “recomputed” (Line 19) 
for the hrst time. Whenever this happens, it will set new the value of Q(x, y, i) equal to 
R(x,y, (1 + y, i)]). Thus, by Lemma 6.3 we will again have 

pShile^x, y, (1 + e)”A’ 2 ^>b+ 2 i^ g) C Q(x, y, i) as demanded by Invariant (12). 

It remains to show that at the last time the algorithm has called Refresh(x^, y', i + l) 
and sets (x',y') as an (i + l)-parent of (x,y), we have 

y, (1 + e)’-A.y++2*^ c y'^ g + gr(+,+.i+l)+2(i+l)^ ^ 

Let vr' denote the shortest path from x' to y' in y', i + 1)] computed at 

the beginning of the last execution of Refresh(x', y', i + l). To enhance the readability of 
this part of the proof we use the abbreviation H = ’^~^^^[Q{x',y',i + 1)]. By the 

if-condition in Algorithm 2 we know that w^ir', H) = dn^x', y') < (1 + eY^^'. 

Consider some node v G y, (1 + e)'’(^>?^'*)+2*^ G) which means that v lies on a path 

vr from x to y in G that has at most Shi/e edges and weight at most (l+e)’'(*’2^’*)+2*. Remember 
that X and y are consecutive i-centers on tt'. Let and tt^ denote the subpaths of vr' from x' 
to X and y to y', respectively. The concatenation vr" = vr'j^ o vr o vr^ is a path from x to y in G. 
We will show that vr" has at most 8/ij+i/e edges and weight at most (1 + 

This then proves that all nodes on vr are contained in R^^*+i/^(x, y, + G). 

As 2hi < hj+i, the number of edges of vr is at most 8 hi/e < 4hj+i/e. Remember that the 
edge weights of H are multiples of e(l + e)”^^ Thus, each edge of H has weight 
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at least e(l + the weight of vr' in H is at most (1 + the 

number of edges of ir' is at most 

(1 + ^Y{x',y',i+l)+2}^.^^ ^ (1 + e)^h,+i ^ 4 
e(l + e)»-(a:'y.*+i) e - e ' 

It follows that the number of edges of vr" is at most 4/ij+i/e + Ahi+i/e = S/ij+i/e. 

By Invariant (II) we have dQ{x,y) > (1 + and thus 

u;G(vr) < (1 + < (l + ef^d^^{x,y) . 

By the initial random sampling of i-centers, every shortest path consisting of hj — 1 edges 
contains an i-center whp (Lemma 2.11). Thus, the subpath of vr' from x io y has at most hi 
edges. Since tt' is a shortest path in H we have dH{x,y) = d^{x,y). It follows that 

wciT^) < (1 + e)^*dG‘(x, y) < (1 + e)^*c^G[Q(x'y,i+i)](^’ y) < (1 + e)^'d^H{x, y) 

= (1 + ef"dH{x,y). 

Furthermore we have duix',y') = dnix',x) + dH{x,y) + dn (and in particular dH{x,y) < 
dnix' ,y')). As wg{tti) < duix' ,x) and wg{t^ 2 ) < dH{y-,y') the weight of the path vr" in G 
can be bounded as follows: 

^^^G(7^") = wciT^'i) + WG{'r^2) + wg{t^) 

< dnix', x) + dniy, y') + wg{t^) 

< dnix, y) - dnix, y) + wg{t^) 

< dnix, y') - dnix, y) + (1 + ef'-dnix, y) 

= dnix', y') + ((1 + e)^* - l)dH{x, y) 

< dnix, y') + ((1 + e)^* - l)dH{x , y) 

= (1 + ef^-duix',y') 

< (1 + ef\l + e)K^'V.*+i)+2 

_ ^yix',y',i+l)+2{i+l) _ Q 

Lemma 6.5. For every 1 <i <k and every active pair of i-centers {x,y) G Ai we have 
dG{x, y) < (1 + whp after the algorithm has finished its updates. 

Proof. The proof is by induction on i. If l(x,y,i) < \Bi\, then there is some Lhub z that 
links X to y. By the triangle inequality we then have 

dG{x, y) < dG{x, z) + dG{z, y) 

< dQSIn/e,r(x,y,i){x, Z) + d^Shi/e,r(x ,y ,i) {z, y) 

< (X _|_ yr(x,y,i)+2i+l _ 

If l{x, y, i) = \Bi\ + 1, then consider the last time the algorithm has called Refresh(x, y, i). 
Let G' and Q'{x,y,i) denote the versions of G and Q{x,y,i) at the beginning of the refresh 
operation and let vr be the shortest path from x to y in {G')^^xi^^yA^Q'(^x,y,i)] computed 
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by the algorithm. To enhance readability we set H = y, i)] in this proof. 

The algorithm ensures that dH{x,y) < (1 + as otherwise, by the if-condition in 

Line 7, {x,y) would have either become inactive (i.e., removed from Ai) or the algorithm 
would have called Refresh(x, y, i) again. Let vi,... ,vi denote the {i — l)-centers in order 
of appearance on tt, i.e., dH{x,y) = J2i<j<i-i dH{vj,Vj+i). 

If i = 1, then remember that the set of 0-centers is equal to V. Note that no edge {u,v) 
of TT has been deleted from G since the last time the algorithm has called Refresh( x, y, i) 
because (x, y) is a 1-parent of {u,v). Therefore all edges of tt still exist in G. Since r{x,y,i) 
is non-decreasing this means that dcix, y) < dH{x,y) < (1 -I- as desired. 

If i > 2, then let 1 < j < / — 1. First of all, note that the value of r{vj,Vj^i,i) has not 
changed since the last time the algorithm has called Refresh(x, y, i ) because otherwise 
after Line 8, the algorithm would call Refresh(x, y, i ) as (x,y) is an i-parent of { vj , Vj ^ i ). 
Therefore, by Invariant (II) of Lemma 6.4, we have d^r^ {vj,Vj^i) > (1 -|- 
Note that H has the same nodes and edges as G' and the weight of each edge in H is at 
least its weight in G'. Therefore dj}~^{vj,Vj+i) > Xj+i) and thus dff~^{vj,Vj+i) > 

(1 -|- gy initial random sampling of {i — l)-centers, every shortest path 

consisting of hi-i — 1 edges contains an (i — l)-center whp (Lemma 2.11). Thus, the subpath of 
TT from Vj to Vj+i has at most hi-i edges which means that duivjiVj+i) = d^~^{vj^Vj+i) > 
(1 _|_ 

By the induction hypothesis we have dcivj^Vjj^i) < (1 -|- in 

the current version of G, which by the argument above implies that dcivj^Vj+i) <(1-|- 
e)^{*-^)'i'^d/f(uj,Uj_|_i). Thus, by the triangle inequality we get 

dcix^y) < da{vj,Vj+i) 

= {1 + eY"~^dH{x,y) 

< (l + e ) 2 *- i(l + e )^(^’^’*)+2 

= {1 + e)2*+^ 


as desired. □ 

Lemma 6.6. After the algorithm has finished its updates we have dcisfi) < S{s,t) < 
{1 + eY~^^dG{s,t) whp. 

Proof. Consider the last time the algorithm has called Refresh(s, t, /c -|- 1). Let G' denote 
the versions of G at the beginning of the refresh operation and let tt be the shortest path from 
X to y in G computed by the algorithm. Note that 6{s, t) = (1 -|- eY^~^^dG’{s, t). As distances 
are non-decreasing under deletions in G we trivially have S{s,t) < (1 -|- eY^^dG{s,t). 

Let vi,...,vi denote the fc-centers in order of appearance on tt , which means that 
dG'{x, y) = J2i<j<i-i dG'{xj,Vj+i), and let 1 < j < / — 1. The value of r{vj,Vj+i, k) has not 
changed since the last time the algorithm has called Refresh(x, y, i) because otherwise after 
Line 8, the algorithm would call Refresh(s, t, k) as (s,t) is a (A: -|- l)-parent of {vj,Vj+i). 
By Lemma 6.4, we have d^,{vj,Vj^i) > (1 -|- and by Lemma 6.5 we have 

dG{vj,Vj+i) < (1 -|- e)^K>^i+i>^)+2^+i, g fgiiQ^g g^^at dG{vj,Vj+i) < (1 -|- eY^~^^dGYvj,Vj+i). 
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By the initial random sampling of /c-centers, every shortest path consisting of /ifc — 1 
edges contains a /c-center whp (Lemma 2.11). Thus, the subpath of vr from vj to vj+i has 
at most hk edges which means that dG'{vj,Vj+i) = Uj+i). Thus, by the triangle 

inequality we get 


dG{s,t) < dG{vj,Vj+i) < {1 + ef^+^dG'{vj,Vj+i) 

= {I + €)‘^’"^^dG'{s,t) = 6{s,t). □ 


Lemma 6.7. For all 0 < x < 1 and all y > 0, 


Proof. Let e denote Euler’s constant. We will use the following well-known inequalities: 
(1 -f l/zY < e (for all z > 0), < 1/(1 — z) (for all 2 ; < 1), and 1/(1 — z) < 1 + 2z (for all 

0 < z < 1/2). We then get: 


< 62 < 


1 - - 
2 


<1 + X . 


□ 


6.4 Running Time 


The running time analysis follows similar arguments as for the s-t reachability algorithm 
in Lemma 5.9. For every pair of i-centers (x,y), the algorithm never decreases r{x,y,i) 
and {x,y) causes a refresh only if r{x,y,i) increases, which happens OilogW/e) times. 
Furthermore it can be argued in a straightforward way that for every 1 < i < k+1 and every 
pair of i-centers {x,y) the graph G[Q{x,y,i)] has at most mm.{mlbi,n^ jhl) edges whp. Now 
the only differences compared to the running time analysis of the s-t reachability algorithm 
are the number of path unions of active pairs of centers each node is contained in and the 
time needed for maintaining the hub links, both of which are analyzed below. 


Lemma 6.8. Let 1 <i <k and consider a pair of active {i -|- l)-centers {x',y') and their 
i-children ixj,yj)i<j<i (which are active i-centers). Then, for every node v, there are at 
most q = 2^*+^ |'logi_,_g(niy)] pairs of i-children {xj,yj) of {x',y') such that v G Q{xj,yj,i). 

Proof. We show that at the last time the algorithm has called Refresh(x/ y', i + 1) (where 
it determined the current i-children of {x',y')) there are at most q = 2^*+^ |'log]^_,_£(nlT)] 
pairs of i-children {xj,yj) of {x',y') such that 


V G V{Xj,yj, (1 + ^ ^ 


Note that, for each 1 < i < /, r{xj,yj,i) has not changed since the last refresh (as such 
a change implies refreshing (x',y')). Thus, for each 1 < i < i, Q{xj,yj,i) is a subset of 
the path union above, which means that the lemma will be implied by this claim. Now 
we actually prove the following slightly stronger claim: for every every node v and every 
0 < r < [log;^^^(nVF)J , there are at most q' = 2^*+^ pairs of i-children {xj,yj) of {x',y') 
such that 


u G V{xj,yj, (1 + + 1)]) and 4‘(x,-,%) > (1 + ef 
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In this proof we use the abbreviation H = ,y',i + 1)]. Suppose that v is 

contained in q' > path unions V{xj,yj,{l + , H) of i-children {xj,yj) such 

that dQ{xj,yj) > (1 + e)'". Let ji, ... ,jq> be the corresponding indices and assume without 
loss of generality that ji < j 2 ■ ■ ■ < jq'- By the initial random sampling of i-centers, every 
shortest path consisting of hj — 1 edges contains an z-center whp (Lemma 2.11). Therefore, 
for every 1 < j < I, there is a shortest path between xj and yj in H with at most hi 
edges, i.e., dH{xj,yj) = d^{xj,yj). Furthermore, for every j G {ji, ■ ■ ■, jq'}^ have 
dQ{xj,yj) > (1 + e)'" by our assumption and thus 

dH{xj^,yj^,) > dH{xj,yj)= d’]^{xj,yj) 

> Y. E (i+er = 9'(i+«)’' (10) 

We now derive an upper bound on dnixj ^, yj^,) contradicting this lower bound. Since v is 
contained in the path union V{xj, yj, (1 + e)’’+2*+i^ //)^ know by the definition of the path 
union that v lies on a shortest path from xj^ to yj^ in H of weight at most (1 + e)'’+2*+i and 
thus dn^Xj^^v) < (1 + e)^+2*+i. xhe same argument shows that dH{v,yj^,) < (1 + e)^+2*+i. 
By the triangle inequality we therefore have 

dH{xj,,yj^,) < dH{xj,,v) + dH{v,yj^,) < 2(1 + e)’'+2*+i (n) 

By combining Inequalities (10) and (11) we get q' < 2(1 + which is a contradictory 

statement for q' > 2^*+^ because e < 1. Thus, v is contained in at most q' = 2^*"’“^ path 
unions V{xj,yj, (1 + e)'’, H) of i-children {xj,yj) such that dQ{xj,yj) > (1 + e)'’. □ 

We now bonnd the time needed for maintaining the hnb links as follows. For every z-hnb 
z (with 1 < i < k) and every 0 < r < logi_^_^{nW) we maintain both an incoming and an 
outgoing ES-tree up to depth (1 + g)^+2*+i In every edge weight is a 

multiple of p = e^(l + e)^/(8/ij) and thus we can scale down the edge weights by the factor 
1/p and maintain the ES-tree in the resulting integer-weighted graph up to depth 

+ ^ 2 ‘^^+\l + eY ^ 2 ^^+Yl + eY8hi ^ 

p ~ p e‘^{l + eY 

Thus, maintaining these two trees takes time 0{2‘^^mhi/e^) which is 0{2‘^^mn/[cie^)) as 
hi = n/ci. As we have 0(logi_|_^(nIF)) = 0(logVF/e) such trees for every i-hub and there 
are 0{hi) i-hubs in, maintaining all these trees takes time 0(X)i<i<fc 2^^6jmnlog IF/(cie^)). 

We now bound the time needed for maintaining the index l{x,y,i) for every pair of 
z-centers (x, y). First, observe that l{x, y, i) assumes integer values from 1 to \Bi\ (the number 
of z-hubs) and it only decreases (to 0) if r(x, y, i) increases. The index r(x, y, i) on the other 
hand is non-decreasing and assnmes integer valnes from 0 to [log^_,_g(nIF)] = 0{\ogW/ e). As 
\Bi\ = 0{bi), the valne of l{x, y, i) therefore changes 0{bi log W/e) times. As there are 0(2^Cj) 
z-centers, the indices of all pairs of centers together change at most 0(X)i<j<fc 2‘^^bicf log IF/e) 
times. It remains to bonnd the time spent in total for calls of the form UpdateHubLinks(x, 
y, i, r) in which the index /(x, y, i) does not increase bnt the algorithm still spends constant 
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time for checking whether l{x, y, i) should increase. In such a case we charge the running time 
to the pair of i-centers {x,y). Note that the call UpdateHubLinks(x, y, i, r) has either 
happened because (x, y) is made an i-child of some pair of z + 1-centers or because (x, y) is an 
active pair of z-centers (i.e., in Ai) after the deletion of a node. We only have to focus on the 
second case because in the first case the charge of 0(1) on (x, y) can be neglected. As argued 
above, at any time, there are at most 0{q^ \ogij^^{nW)) active z-centers. Furthermore, there 
are at most m deletions in G. Therefore the total time needed for maintaining all indices 
/(x, y, z) is 0(X)i<j<fe 2‘^^bicf log Wje -h y^mlog W + e). 

Putting everything together, the total running time of our algorithm using k layers is 


bicf log W ^ bimn log W 2 k . I m n 


o E E 








+ q m- mm ^, vr 


+ E 


Ca mm 






i+i. 


where q = 2^^+^ [log^^^ W]. Assuming that W < and e > l/log'^ n and setting 

k = [log^/^nj we get 2^q = We now simply set the 

parameters bi <n and Ci < n for each 1 < z < A: in the same way as in Lemma 5.9 to obtain 
the same asymptotic running time (in terms of polynomial factors) as for the s-t reachability 
algorithm. 

Theorem 6.9. For every VF > 1 and every 0 < e < 1 such that W < and e > 1/log'^rz 

for some constant c, there is a decremental (1 -|- e)-approximate stSP algorithm with constant 
query time and total update time 

0(min(m5/^nV2+«(i),„i2/3^4/3+o(i)^^ ^ 0 {mn^n+oW) 


that is correct with high probability against an oblivious adversary. 

Similar to the s-t reachability algorithm we can extend the (1 -|- e)-approximate stSP 
algorithm in the following ways. 

Corollary 6.10. For every IT > 1 and every 0 < e < 1 such that W < 2*°®'""' and 
€ > l/log'^rz for some constant c, there is a decremental (1 -|- e)-approximate algorithm for 
maintaining shortest paths between p source-sink pairs with constant query time and total 

iiT)flnfp fiTinp 

0(min(pV2^5/4^1/2+o(l)^^l/3^2/3^4/3+o(l)^^ 

that is correct with high probability against an oblivious adversary. 

Corollary 6.11. For every IT > 1 and every 0 < e < 1 such that W < 2*°®'"” and 
€ > l/log'^ n for some constant c, there is a decremental (1 -|- e)-approximate SSSP algorithm 
with constant query time and total update time 

0(min(m^/®n2/3+«(P,m3/4n5/^+'’(^))) = 


that is correct with high probability against an oblivious adversary. 
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7 Faster Single-Source Reachability in Dense Graphs 


In this section we first introduce a path union data structure that is more efficient than the 
naive approach for repeatedly computing path unions between one fixed node and other 
variable nodes. We then show how to combine it with a multi-layer approach to obtain a 
faster decremental single-source reachability algorithm for dense graphs. 

7.1 Approximate Path Union Data Structure 

In the following we present a data structure for a graph G undergoing edge deletions, a 
fixed node x, and a parameter h. Given a node y, it computes an “approximation” of the 
path union 'P(x, y, h, G). Using a simple static algorithm the path union can be computed 
in time 0{m) for each pair {x,y). We give an (almost) output-sensitive data structure for 
this problem, i.e., using our data structure the time will be proportional to the size of the 
approximate path union which might be o{m). Additionally, we have to pay a global cost of 
0 {m) that is amortized over all approximate path union computations for the node x and all 
nodes y. This will be useful because in our reachability algorithm we can use probabilistic 
arguments to bound the size of the approximate path unions. 

Proposition 7.1. There is a data structure that, given a graph G undergoing edge deletions, 
a fixed node x, and a parameter h, provides a procedure ApproximatePathUnion such 
that, given sequence of nodes yi,... ,yk, this procedure computes sets Fi,... Fk guaranteeing 
V{x,y,h,G) F Fi C V{x,y, (logm -|- 3)h,G) for all 1 < i < k. The total running time is 
o{E l<i<k \^i \ + "^)- 

7.1.1 Algorithm Description 

Internally, the data structure maintains a set R{x) of nodes, initialized with R{x) = V, such 
that the following invariant is fulfilled at any time: all nodes that can be reached from x by 
a path of length at most h are contained in R{x) (but R{x) might contain other nodes as 
well). Observe that thus R{x) contains the path union V{x,y,h,G) for every node y. 

To gain some intuition for our approach consider the following way of computing an 
approximation of the path union V{x, y, h, G) for some node y. First, compute Bi = {v G 
R{x) I y) ^ using a backward breadth-first search (BFS) to y in G[R{x)], the 

subgraph of G induced by R{x). Second, compute F = {v G R{x) \ dQ[Bi]ix,v) < h} using 
a forward BFS from x in G[Bi]. It can be shown that T’{x, y,h,G) C. F C V{x, y, 2h, G).^^ 
Given Bi, we could charge the time for computing F to the set F itself, but we do not know 
how to pay for computing Bi as Bi\F might be much larger than F. 

Our idea is to additionally identify a set of nodes X F {v gV \ dcix, v) > h} and remove 
it from R(x). Gonsider a second approach where we first compute Bi as above and then 
compute B 2 = {v G R{x) \ dG[R(x)](u,y) < 2/i} and F = {v G R{x) \ dQ[B^^{x,v) < h}. It 
can be shown that V{x, y,h,G) F F GL V{x, y, 3h, G). Additionally, all nodes in A = Bi \ F 
are at distance more than h from x and therefore we can remove X from R{x). Thus, we 

'^^Indeed, F might contain some node v with dG{x,v) = h and dG{v,y) = h, but it will not contain any 
node w with either dG(x,w) > h or dG(w,y) > h. 
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can charge the work for computing Bi and F io X and F, respectively.^^ However, we now 
have a similar problem as before as we do not know whom to charge for computing i? 2 . 

We resolve this issue by simply computing = {u E R{x) \ 2/) < 

increasing values of i until we arrive at some i* such that the size of Bi* is at most double 
the size of We then return F = {u E R{x) \ dQ[Bi]{x,v) < h} and charge the time 

for computing Bi to X = Bi^i \ F and F, respectively. As the size of Bi can double at 
most O(logn) times we have V{x,y,h,G) C F C V{x,y,0{hlogn),G), as we show below. 
Procedure 3 shows the pseudocode of this algorithm. Note that in the special case that x 
cannot reach y the algorithm returns the empty set. In the analysis below, let i* denote the 
final value of i before Procedure 3 terminates. 


Procedure 3: ApproximatePathUnion(7/) 


// All calls of ApproximatePathUnion(?/) use fixed x and h. 

1 Compute Bi = {v G R{x) \ dQ[ji{x)]ix, y) < h} // backward BFS to y in subgraph 

induced by R(x) 

2 for i = 2 to [logm] + 1 do 

Compute Bi = {v G R{x) \ dG[R(x)](f, y) ^ // backward BFS to y in 

subgraph induced by R(x) 
if \E[Bi]\ <2|F[F,_i]| then 

Compute F = {v G Bi \ dQ[Bi]{x,v) < h} // forward BFS from x in 
subgraph induced by Bi 
X G- Bi_i \ F, R{x) G- R{x) \ X 

return F 


7.1.2 Correctness 

We first prove Invariant (I): the set R{x) always contains all nodes that are at distance at 
most h from x in G. This is true initially as we initialize R{x) to be V and we now show 
that it continues to hold because we only remove nodes at distance more than h from x. 

Lemma 7.2. If R{x) C {u E P | dG{x,v) < h}, then for every node v G X removed from 
R{x), we have dc^x^v) > h. 

Proof. Let v G X = Fj*_i \ F and assume by contradiction that dcixjv) < h. Since 
v G Fj*_i we have dG[R(a:)](^) 2/) ^ “ 1)^- Now consider the shortest path vr from x to u 

in G, which has length at most h. By the assumption, every node on tt is contained in 
G[R{x)]. Therefore, for every node v' on tt, we have dG[R(x)]W,v) < h and thus 

dG[R{x)]{v',y) < dG[R{x)]{v\v) + dG[R{x)]{v,y) <h + {i* - l)h < i*h 

which implies that v' G Bi*. Thus, every node on vr is contained in Bi*. As tt is a path 
from X to r of length at most h it follows that dG[B^*]ix, v) < h. Therefore v G F, which 
contradicts the assumption v G X. □ 

'^^Note that in our first approach removing Bi\ F would not have been correct as F was computed w.r.t 
to G[Bi\ and not w.r.t. G[B 2 \. 
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We now complete the correctness proof by showing that the set of nodes retnrned by the 
algorithm approximates the path union. 

Lemma 7.3. Procedure 3 returns a set of nodes F such that V{x^ y,h,G)CFC V{x, y, (log m+ 
3)h,G). 

Proof. We first argue that the algorithm actually returns some set of nodes F. Note that in 
Algorithm 3 of the algorithm we always have as i?j_i C Bi. As E[Bi] is 

a set of edges and the total nnmber of edges is at most m, the condition 
therefore must eventually be fulfilled for some 2 <i < [logm] + 1. 

We now show that V{x, y, h, G) C F. Let v G V{x, y, h, G), which implies that v lies on 
a path TT from x to y of length at most h. For every node v' on tt we have dcix^v') < h, 
which by Invariant (I) implies v' G R{x). Thus, the whole path vr is contained in G[i?(x)]. 
Therefore 2/) ^ h for every node v' on vr which implies that tt is contained in 

G[Bi*]. Then clearly we also have v) <h which implies v G F. 

Finally we show that F C V{x, y, (log m + 3)/i, G) by proving that dcix, v) + dciv, y) < 
(logm + 3)h for every node v G F. As G[Bi*] is a subgraph of G, we have dG{x,v) < 
dG[Bi*]{x,v) and dG{v,y) < dG[Bi*]{v, y)■ By the definition of F we have dG[Bi*]{x,v) < h. 
As F G Bi* we also have dG[B^*]{v,y) < i*h < ([logm] + l)h < (logm + 2)/i. It follows that 
dG{x, v) + dG{v, y) < h + (log m + 2)h = (log m + 3)h. □ 

7.1.3 Running Time Analysis 

To bonnd the total rnnning time we prove that each call of Procedure 3 takes time proportional 
to the number of edges in the returned approximation of the path union plus the number of 
edges incident to the nodes removed from R{x). As each node is removed from R{x) at most 
once, the time spent on all calls of Procedure 3 is then 0{m) plus the sizes of the subgraphs 
induced by the approximate path unions returned in each call. 

Lemma 7.4. The running time of Procedure 3 is 0{\E[F]\ + \E[X, R{x)]\ + \E[R{x),X]\) 
where F is the set of nodes returned by the algorithm, and X is the set of nodes the algorithm 
removes from R{x). 

Proof The running time in iteration 2 < j < i* — 1 is 0{\E[Bj]\) as this is the cost of 
the breadth-first-search performed to compute Bj. In the last iteration i*, the algorithm 
additionally has to compute F and X and remove X from R{x). As F is computed by a 
BFS in G[Bi*] and X C Bi*-i C Bi*, these steps take time 0{\E[Bi*]\). Thus the total 
running time is 0{Y,i<j<i* \E[Bj]\). 

By checking the size bound in Line 4 of Procedure 3 we have |Fl[i3j]| > 2|Fl[i?j_i]| for 
all 1 < j < P — 1 and |Fl[i3i*]| < 2\E[Bi*-i]\. By repeatedly applying the first inequality it 
follows that < 2\E[Bi*-i]\. Therefore we get 

^ \E[B,]\= \E[B,]\ + \E[B,*]\<2\E[B,*.,]\+2\E[B,*.i]\=4\E[B,*.i]\ 

and thus the running time is 0{\E[Bi*-i]\). Now observe that by A = \F we have 

Bi*-i G X U F and thus 

E[Bi*_i] G E[F] U E[X] U E[X, F] U E[F, X] G E[F] U E[X, R{x)] U E[R{x),X] . 
Therefore the running time is 0(|F1[F]| -|- \E[X, R{x)]\ -|- \E[R{x), X]\). □ 
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7.2 Reachability via Center Graph 

We now show how to combine the approximate path union data structure with a hierarchical 
approach to get an improved decremental reachability algorithm for dense graphs. The 
algorithm has a parameter 1 < A: < log re and for each 1 < i < k a, parameter Ci < n. We 
determine suitable choices of these parameters later in this section. For each 1 < i < k — 1, our 
choice will satisfy q > Cj+i and Cj = 0(ci+i). Furthermore, we set hi = (3+logrre)*“^re/ci for 
1 <i <k. At the initialization, the algorithm determines sets of nodes Ci A (72 A • • • D 

such that s, t G (7i as follows. For each 1 < i < A:, we sample each node of the graph with 

probability aci In re/re (for a large enough constant a), where the value of q will be determined 
later. The set Ci then consists of the sampled nodes, and if i < /c — 1, it additionally contains 
the nodes in (7j+i. For every 1 < i < /c we call the nodes in Ci i-centers. In the following we 
describe an algorithm for maintaining pairwise reachability between all 1-centers. 

7.2.1 Algorithm Description 

Data Structures. The algorithm uses the following data structures: 

• For every i-center x and every i < j < k an approximate path union data structure 
(see Proposition 7.1) with parameter hj. 

• For every A:-center x an incoming and an outgoing ES-tree of depth in C. 

• For every pair of an i-center x and a j-center y such that I := max(i,j) < A: — 1, a set 

of nodes Q{x, y, 1) C V. Initially, Q{x, y, 1) is empty and at some point the algorithm 
might compute Q{x, y, 1) using the approximate path union data structure of x. 

• For every pair of an i-center x and a j-center y such that I := max(i, j) < A; — 1 an 
ES-tree of depth hi from x in Q{x,y,l). 

• Eor every pair of an i-center x and a j-center y such that I := max(i, j) < A: — 1 a set 
of (/ + l)-centers certifying that x can reach y. 

Certified Reachability Betweeu Ceuters (Liuks). The algorithm maintains the fol¬ 
lowing limited path information between centers, called links, in a top-down fashion. Let 
X be a A:-center and let y be an i-center for some 1 < i < A: — 1. The algorithm links x 
to y if and only if y is contained in the outgoing ES-tree of depth hk of x. Similarly the 
algorithm links y to x if and only if y is contained in the incoming ES-tree of depth hk of 
X. Let X be an i-center and let y be a j-center such that I := max(i, j) < A: — 1. If there is 
an {I + l)-center 2 ; such that x is linked to 2 and z is linked to y, the algorithm links x to 
y (we also say that z links x to y). Otherwise, the algorithm computes Q{x,y,l) using the 
approximate path union data structure of x and starts to maintain an ES-tree from x up 
to depth hi in G[Q{x, y, 1)]. It links x to y if and only if y is contained in the ES-tree of x. 
Using a list of centers z certifying that x can reach y, maintaining the links between centers 
is straightforward. 
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Center Graph. The algorithm maintains a graph called center graph. Its nodes are the 
1-centers and it contains the edge (x, y) if and only if x is linked to y. The algorithm 
maintains the transitive closure of the center graph. A query asking whether a center y is 
reachable from a center x in G is answered by checking the reachability in the center graph. 
As s and t are 1-centers this answers s-t reachability queries. 

7.2.2 Correctness 

For the algorithm to be correct we have to show that there is a path from s to t in the center 
graph if and only if there is a path from s to t in G. We will in fact show in more generality 
that this is the case for any pair of 1-centers. 

Lemma 7.5. For every pair of an i-center x and a j-center y, if x is linked to y, then there 
is a path from x to y in G. 

Proof. The proof is by induction on I = max(z,j). As x is linked to y, one of the following 
three cases applies: 

1. j = k and x is contained in the incoming ES-tree of depth hk of y in G. 

2 . i = k and y is contained in the outgoing ES-tree of depth hk of x in G. 

3. There is an {I -\- l)-center z such that x is linked to z and z is linked to y. 

4. y is contained in the ES-tree of depth hi of x in G[Q(x, y, /)]. 

In the first two cases there obviously is a path from x to y in G by the correctness of the 
ES-tree. In the third case we may apply the induction hypothesis and find a path from x to 
z as well as a path from z to y in G. The concatenation of these paths is a path from x to y 
in G. In the fourth case we know by the correctness of the ES-tree that there is a path from 
X to y in G[Q(x, y, 1)] and therefore also in G. □ 

Lemma 7.6. For every pair of an i-center x and a j-center y, if dG{x,y) < hi, then x is 
linked to y. 

Proof. Set I = max(i,y). If I = k, then assume that I = i (the proof for I = j is symmetric). 
Since dQ{x,y) < hk, y is contained in the outgoing ES-tree of depth hk of x by the correctness 
of the ES-tree. Thus, x is linked to y. 

If / < fe — 1 and there is an [I -\- l)-center z such that x is linked to z and z is linked 
to y, then also x is linked to y. If this is not the case, then the algorithm has computed 
Q{x,y,l) using the approximate path union data structure and at that time we have 
V{x,y,hi,G) C Q(x,y,l) by Proposition 7.1. By Lemma 2.8 the set Q{x,y,l) contains 
V{x, y, hi, G), also in the current version of G. As dcix, y) < hi, all nodes on the shortest 
path from x to y in G are contained in V{x, y, hi, G) and thus in Q{x, y, 1). Therefore the 
ES-tree from x up to depth hi in G[Q{x,y,l)] contains y which means that x is linked to 
y. □ 

Lemma 7.7. For every pair of 1-centers x and y, there is a path from x to y in the center 
graph if and only if there is a path from x to y in G. 
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Proof. Assume that there is a path from x to y in the center graph. Every edge {x',y') in 
the center graph (where x' and y' are 1-centers) can only exist if the algorithm has linked x' 
to y'. By Lemma 7.5 this implies that there is a path from x' to y' in G. By concatenating 
all these paths we obtain a path from x to y in G. 

Now assume that there is a path from x to y in G. Consider the shortest path vr from x to 
y in G and any two consecutive 1-centers x' and y' on tt. By the initial random sampling of 
1-centers, every shortest path consisting of hi — 1 edges contains a 1-center whp (Lemma 2.11) 
and thus dcix', y') < hi. By Lemma 7.6 the algorithm has linked x' to y' and thus there is 
an edge from x' to y' in the center graph. As snch an edge exists for all consecutive 1-centers 
on TT, there is a path from x to y in the center graph. □ 

7.2.3 Running Time Analysis 

The key to the efficiency of the algorithm is to bound the size of the graphs Q(x, y, 1). 

Lemma 7.8. Let x he an i-center and let y be a j-center such that I := max(i, j) < k — 1. 
If X is not linked to y by an {I l)-center, then Q{x, y, 1) contains at most n/q+i nodes with 
high probability. 

Proof. Suppose that Q{x,y,l) contains more than n/q+i nodes. Then by the random 
sampling of centers, Q{x, y, 1) contains an (Z-l-l)-center z with high probability by Lemma 2.11. 
By Proposition 7.1 we have dc{x, z)-\-da{z,y) < (3-|-log m)/i; = and thus dG(x, z) < 
and dciz, y) < It follows that x is linked to z and z is linked to y by Lemma 7.6. But 
then X is linked to y, contradicting our assumption. □ 

With the help of this lemma we first analyze the running time of each part of the 
algorithm and argue that our choice of parameters gives the desired total update time. 


Parameter Choice. We carry out the running time analysis with regard to two parameters 
1 < 5 < c < n which we will set at the end of the analysis. We set 


k 


log(c/fe) 

\/log n • log log n 


+ 1 ) 


Ck = b and c* = = 0(cj+i) for 1 < i < fc — 1. Note that the number of 

i-centers is 0(cj) whp. Observe that 


(3 + logm)^-^ = 0((logn)^) < 0((logn)Vi°g^/i°gi°g’^) 

— ( 9 (' 2 v'logn-log logn^ _ lagn/ logn^ _ ^ 

Furthermore we have 

Cl = (2^^°Sn-loglogn^^-l ^ 2^og{c/b)^ = . b = C 

and by setting k' = (log (c/6))/(\/log n • log log n) we have k <k' + 2 and thus 


^ fi-loglogn^'^ _ 2 \/logn-loglog_ 

Remember that hi = {2> + logm)*“^n/ci for 1 <i <k. Therefore we have hi = 0(n/ci) = 
0{n/c). 


k'+l 
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Maintaining ES-Trees. For every /c-center we maintain an incoming and an outgoing 
ES-tree of depth hk, which takes time 0{mhk). As there are 0{ck) /c-centers, maintaining 
all these trees takes time 0 {ckmhk) = 0 {bmn/c). 

For every i-center x and every j-center y such that I := max{i,j) < k — 1, we maintain 
an ES-tree up to depth hi in G'[Q(x, y, /)]. By Lemma 7.8 Q{x, y, 1) has at most u/q+i nodes 
and thus G[(5(a:, y,/)] has at most edges. Maintaining this ES-tree therefore takes 

time 0{{'n?■ hi) = 0{n‘^/cf_^_^{n/ci)) = 0{n^/{cicfj_i)). In total, maintaining all these 
trees takes time 


O 


ciCj 

l<i<k-ll<j<i 


ClC?+i j 


6 


E E 

l<i<k-l l<j<i 


CiCin^ 

Ci+lCiCk 


= 6 


E E 

l<i<fe-l l<i<i 




Computing Approximate Path Unions. For every i-center x and every i < j < k 
we maintain an approximate path union data structure with parameter hj. By Proposi¬ 
tion 7.1 this data structures has a total running time of 0{m) and an additional cost of 
0{\E[Q{x,y,j\)\) each time the approximate path union Q{x,y,j) is computed for some 
j-center y. By Lemma 7.8 the number of nodes of Q{x,y,j) is njcj+i with high probability 
and thus its number of edges is Therefore, computing all approximate path unions 

takes time 


O 


E E 

vl<i<k—l i<j<k 


n 


CiTfl H- CiCj 




= o{ E E Cl"" + 


= ^ E E Cl"' + 


.1 i<j<k 


Cin 

Ck 


yl<i<k—l i<j<k 
2 


ciCj-n^ \ 

j 


= 0{k^cim -|- k"^cin^/ck) = 0{cm -|- crc jh) . 


Maintaining Links Between Centers. For each pair of an i-center x and a y-center y 
there are at most 0{ci^i) (/-|-l)-centers that can possibly link x to y. Each such (/-|-l)-center 
is added to and removed from the list of (/ -|- l)-centers linking x to y at most once. Thus, 
the total time needed for maintaining all these links is 0 {J2i<i<k-iJ2i<j<i'^i(^jCi+i) = 
0{k‘^c\) = O(c^). 


Maintaining Transitive Closure in Center Graph. The center graph has 0{ci) nodes 
and thus 0(cf) edges. During the algorithm edges are only deleted from the center graph and 
never inserted. Thus we can use known 0{mn)-time decremental algorithms for maintaining 
the transitive closure [RZ08, Lacl3] in the center graph in time 0(cf) = O(c^). 

Total Running Time. Since the term cr? jh is dominated by the term n^/6, we obtain 
a total running time of O (^bmnjc + r? jb -|- cm -|- (?). By setting b = jm?/'^ and c = 

running time is jm) and by setting b = and 

c = the running time is 0 {jr?l’^ 
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7.2.4 Decremental Single-Source Reachability 

The algorithm above works for a set of randomly chosen centers. Note that the algorithm 
stays correct if we add any number of nodes to Ci, thus increasing the number of 1-centers 
for which the algorithm maintains pairwise reachability. If the number of additional centers 
does not exceed the number of randomly chosen centers, then the same running time bounds 
still apply. Thus, from the analysis above, we obtain the following result. 

Theorem 7.9. Let S CV be a set of nodes. If\S\ < , then there is a decremental 

algorithm for maintaining pairwise reachability between all nodes in S with constant query 
time and a total update time of + n^/m) that is correct with high probability 

against an oblivious adversary. If \S\ < , then there is a decremental algorithm 

for maintaining pairwise reachability between all nodes in S with constant query time and 
total update time that is correct with high probability against an 

oblivious adversary. 

Using the reduction of Theorem 4.1 this also gives us a single-source reachability algorithm. 

Corollary 7.10. There are decremental SSR and SCC algorithms with constant query time 
and total update time 

0(^2/3„4/3^ ^3/7^12/7) 

that are correct with high probability against an oblivious adversary. The total update time 
of the SCC algorithm is only in expectation. 

Corollary 7.11. There is a decremental SSR algorithm with constant query time and total 
update time 

+ m3/7„i2/7) 

that is correct with high probability against an oblivious adversary. 

Proof. First, set ci = and observe that by Theorem 7.9 we have an algorithm 

that can maintain reachability from a source to c\ sinks with a total update time of 
+n‘^/m). By Theorem 4.1 this implies a decremental single-source reachability 
algorithm with a total update time of 0{nn?/^n^/^+n^/m+mn/ci). The same argument gives 
a decremental single-source reachability algorithm with a total update time of + 

P mnfc 2 ) where C 2 = m}!'^. 

If m > n®/® we have jm < and mnjci = jnfl'^ < If m < 

we have and mnjc^ = < rrpl'^rtf’ll. Thus, by running the 

first algorithm if m > vPP and the second one if m < vPP we obtain a total update time of 
0{mfPn^P - 1 - (Note that rr?!'^n}’^!'^ < if and only if m > n ®/®). □ 

Furthermore, the reduction of Theorem 4.2 gives a decremental algorithm for maintaining 
strongly connected components. 

Corollary 7.12. There is a decremental SCC algorithm with constant query time and 
expected total update time 

0(^2/3^4/3^ ^3/7^12/7) 

that is correct with high probability against an oblivious adversary. 
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8 Conclusion 


In this paper we have presented decremental algorithms for maintaining approximate single¬ 
source shortest paths (and thus also single-source reachability) with constant query time and 
a total update time of o{mn). This result implies that the thirty-year-old algorithm with 
0{mn) total update time is not the best possible. Since it is hard to believe that the running 
time obtained in this paper is tight, it is an important open problem to search for faster 
(and perhaps simpler) algorithms. Note that decremental approximate SSSP in undirected 
graphs can be solved with almost linear total update time [HKN14b]. Thus, it is interesting 
to see if the same total update time can be obtained for the directed case. Getting this 
result for s-t reachability will already be a major breakthrough. 

Given recent progress in lower bounds for dynamic graph problems [PatlO, AVW14, 
HKN+15], it might also be possible that an almost linear update time is not possible for 
decremental approximate SSSP and SSR in directed graphs. However, it might be challenging 
to prove this as all known constructions give the same lower bounds for both the incremental 
and the decremental version of a problem. As incremental SSR has a total update time of 
0{m) we cannot hope for lower bounds for decremental SSR using these known techniques. 
This only leaves the possibility of finding a lower bound for decremental approximate SSR or 
of developing stronger techniques. It also motivates studying the incremental approximate 
SSSP problem which is intuitively easier than its decremental counterpart, but no as easy as 
incremental SSR. 

Finally, we ask whether it is possible to remove the assumption that the adversary is 
oblivious, i.e., to allow an adaptive adversary that may choose each new update or query 
based on the algorithm’s previous answers. This would immediately be guaranteed by a 
deterministic algorithm. 
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